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Abstract. Electrostatic plasma confinement and magnetic electrostatic plasma 
confinement (MEPC) have been studied for four decades. The multiple potential well 
hypothesis, postulated to explain high neutron yields from Hirsch's colliding beam 
experiment, has been supported by several pieces of evidence, but results were 
inconclusive. Magnetic shielding of the grid was developed to reduce the required beam 
current and to prevent grid overheating. Electrostatic plugging of magnetic cusps 
evolved to a similar configuration. Due to low budgets, early MBPC experiments used 
spindle cusps, which are poor for plasma confinement. Later experiments used multipole 
cusps or a linear set of ring cusps, which have larger volumes of field-free plasma. To 
keep the self-shielding voltage drop A<£ 100 kV, the electron density n a in the anode 

gap should be less than about 10 19 m“ 3 . The central plasma density can be an order of 
magnitude higher. The ATOLL toroidal quadrupole had anomalous electron energy 
transport, but the Jupiter-2M linear set of ring cusps achieved a transport rate about a 
factor of two above the classical rate. With near-classical transport, a power gain ratio 
Q~10 is predicted for a reactor with r p = 3m, £ a = 6T, and applied voltage 
a = 400 kV. Besides producing electricity and synthetic fuels, MEPC reactors could be 
used for heavy ion beams sources and neutron generators. The main issues of concern 
for MEPC reactor development are electron transport, plasma purity and electrode 
alignment and voltage holding. 


List of parameters 

(SI Units are used, unless otherwise noted.) 

a = half-width of the anode gap 

A = area 

A [oss - effective area for plasma loss out of magnetic cusps 
A = magnitude of the magnetic vector potential 
B = magnetic field (magnetic induction) 

J5 a = magnetic field in the anode gap 

B b = [2 /i 0 nk(T e + T)] xri = boundary magnetic field 

B c = peak magnetic field at the coil 

B z = point cusp magnetic field 

c = speed of light 

c s = ion acoustic speed 


t This work was supported in part by the US Department of Energy, Office of Energy Research, contract 
DE- AC07-7 6ID01 570. 
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d = grid wire diameter 

D a = electron diffusion coefficient due to collisions with neutral atoms 

D c i = electron diffusion coefficient due to collisions with ions 

£> max = maximum diffusion coefficient 

jD 0 = total collisional diffusion coefficient 

e = electronic charge 

fi = fraction of incident neutrals which are ionized (instead of charge exchanged) 

/« e = fraction of fusion product alpha energy W a transferred to electrons 

f ai = fraction of fusion product alpha energy W a transferred to ions 

F = geometrical factor having dimensions of reciprocal length 

g = spacing between grid wire centres 

h = width of the anode along the magnetic field direction 

I = electron current from electrodes into the plasma 

/ di = electron diffusion current 

/ e = electron injection current 

k = wavenumber = 2k /A 

k = Boltzmann constant 

k = subscript denoting either electrons or ions 

L = length between ends of a linear cusp device 

L, = distance between adjacent ring cusps 

m = electron mass 

m = number of pairs of multipole conductors = N/2 
M = ion mass 

n = plasma electron density 

« a = uncompensated electron density at centre of anode gap 

n B = Brillouin density 

n k = particle density of species k 

n c = central core electron density 

N c = total number of electrons 

n 0 = central plasma density 

«o = initial plasma density (pulsed electrostatic confinement) 

« ia = ion density in the anode gap 
n t = trapped cold electron density in anode region 
N = number of multipole cusp gaps 

P = (t j ./t ii ) 1/2 

P E = electron heating by the edge electric field that affects the central plasma 

Pi = ion heating by plasma waves 

Prad = radiative power loss 

P w = electron power loss to plasma waves 

P 0 = canonical angular momentum 

Q = (fusion power)/ (input power) 

q = Wp C /a>ce 

r = radius 

r a = point cusp anode radius 

r b = plasma boundary radius 

r 0 = virtual electrode radius (electrostatic confinement and Penning trap) 

r Q = radius of plasma boundary surface in point cusp 

r p = approximate plasma radius at magnetic field boundary 
R = accelerating grid radius (electrostatic confinement) 


Magnetic electrostatic plasma confinement 


1541 


R = magnetic mirror ratio 

R = ring cusp anode radius 

S = plasma surface area 

Si = ionization source 

S m = volume-averaged neutral input rate 

T ek = electron temperature (keV) 

T e — electron temperature 

r e0 = central electron temperature 

u = ion speed 

= ion speed in anode gap 
v = electron velocity 

v Vn = electron density-gradient drift velocity 

V = plasma volume 

w = characteristic half-width of the electron distribution in anode region 

W a = average energy acquired along the slope of the potential well by ions 

= electron energy 
= ion energy 

Wioss = average electron energy loss per incident neutral atom due to ionization 
W x = ion energy perpendicular to the radial direction 
x = r/R 

y = *i/r. 

Vk =z k e<f>klT k 

Y = 2 nRIN 

Zk ~ charge number 

Z k =1 for electrons and 1/2 for ions 

a = a parameter depending on TJ T c 

y = fraction of current intercepted by the grid 

8 = the effective half-width of the cusp gap for plasma loss 

8 = a function of p l5 p 2 

A 4> = self-shielding voltage drop in anode gap 

A r — thickness of plasma edge region 

e = (trapped electron pump-out rate)/ (diffusion rate) 

s = effective open fraction of wire mesh grid 

£ g = geometrical open fraction of wire mesh grid 

e 0 = permittivity of free space 

l - global Brillouin ratio 

In A = Coulomb logarithm 

A = wavelength 

Po = permeability of free space 

v ei = electron-ion momentum-transfer collision frequency 

p a - electron Larmor radius in anode gap 

p b = electron Larmor radius at plasma boundary 

p e = electron Larmor radius 

po = electron Larmor radius in point cusp 

Pi = ion Larmor radius 

p! = inner electron turning radius 

p 2 = outer electron turning radius 

(cru) = DT fusion reaction rate parameter 
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Tj. = electron loss time by diffusion across the magnetic field 
r = pulse length 

T k = characteristic collision time for electrons or ions 
Tcond - time for heat loss by conduction 

T diff = characteristic time for electron diffusion across the magnetic field 

T e - electron lifetime or confinement time 

t e _l = electron perpendicular energy confinement time 

t E || = electron parallel energy confinement time 

T ei =* electron-ion collision time 

T eq = electron-ion equipartition time 

^ = ion confinement time 

t 0 = electron collision time 

T, rap = time for central plasma free electrons to be trapped by edge magnetic field 

<f> A - applied voltage 

4> k - potential barrier for electrons or ions 

4> e - electron potential barrier 

<j > j — ion potential barrier 

<t> p = plasma potential = <£; + A<p 

(o c — cyclotron frequency 

o> ce = electron cyclotron frequency 

<w ci = ion cyclotron frequency 

<w pe - electron plasma frequency" 


1. Introduction 

Magnetic electrostatic plasma confinement (MEPC) uses electrostatic fields to 
enhance plasma confinement by open magnetic fields. The purpose of this review is 
to describe: 

(i) the evolution of MEPC from electrostatic plasma confinement and from open 
magnetic confinement systems; 

(ii) MEPC theory; 

(iii) experimental observations; 

(iv) fusion reactor concepts; 

(v) salient issues. 

MEPC devices have also been called ‘electromagnetic traps’, ‘magnetoelectro- 
static traps’, ‘magnetic inertial-electrostatic confinement’, ‘electrostatically plugged 
cusps’, ‘kinetic-electric-magnetic plasma confinement (KEMP)’, ‘electric-magnetic 
confinement’, ‘Polywell tm ’, and ‘modified Penning trap’. Discussions of electrostatic 
field applications in tandem mirrors, bumpy tori and tokamak helicity injection are 
beyond the scope of this review. 

1.1. Electrostatic plasma confinement 

1.1.1. Virtual electrode formation. Lavrent’ev [1] first proposed electrostatic plasma 
confinement on 22 June 1950, and magnetic field enhancement of electrostatic 
plasma confinement in March 1951. In the early 1960s Farnsworth [2] proposed 
electrostatic plasma confinement for the production of fusion reactions, based upon 
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Figure 1. Concentric spherical electrodes used for electrostatic plasma confinement. 

his experience with spherical multipactor vacuum tubes. Lavrent’ev [3] developed a 
theory of charged particle flow and focusing in plane and spherical geometries. 
Budker [4] calculated the fusion reaction rate that could be attained by trapping ions 
in the negative electrostatic charge of an electron beam. Early electrostatic plasma 
confinement experiments [5] consisted of concentric spherical wire mesh electrodes, 
which accelerated charged particles inwards to produce a virtual electrode inside the 
central sphere (figure 1). Charged particles of the opposite sign would be attracted 
into the virtual electrode. For example, if electrons were injected into the sphere, 
they would produce a virtual cathode [6], which would attract positive ions. In order 
to minimize grid current and heating, it is desirable to have the grid mesh nearly 
transparent to the charged particle flows. 

Deflections of injected particles towards nearby grid wires decrease the effective 
open area of the grid by an amount [7] 

sle g = (1 - l/16TOg /2 ) 2 (1) 

where s is the effective open fraction of the grid, e g = (g — d) 2 /g 2 is the geometrical 
open fraction of the grid, g is the spacing between grid wire centres (square mesh) 
and d is the grid wire diameter. For example, if g = 4 mm and d = 0.2 mm, then 
s g = 0.90 and s = 0.865. The circulating current 7 C is related to the current I 
intercepted by the grid wires by the equation 

IJI = 2e/(l- s 2 ). (2) 

For the example case with e = 0.865, the theoretical value of IJI = 6.9. Thus, for a 
highly open grid the circulating current can be about an order of magnitude greater 
than the grid current. 

Lavrent’ev and co-workers [8] studied ion generation by oscillating electrons and 
ion focusing with spherical grids. Probe measurements inside a 5 cm radius spherical 
cathode indicated a strong focusing of injected ions down to radii < lcm. The 
focused ion current density was 12-25 mA cm -2 at 4 kV in 0.13 Pa of argon [9-11]. 
A passing ion density of 10 I7 m -3: was deduced from current probe data [12]. 

Hirsch [13] injected electrons into a spherical wire mesh anode, and used an 
electron beam probe to measure the potential well depth produced by the virtual 
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cathode. He found that deep, negative potential wells were easily produced and 
maintained. He also observed some oscillations of the discharge at frequencies 
-100 MHz. 

For the case of electron injection to produce a virtual cathode, there are two 
great difficulties: the very high currents required [14] and grid wire overheating. 
Better confinement could ■ be obtained for the case of ion injection into the 
accelerating grid, producing a virtual anode. 

Hockney [15] simulated the formation of a virtual electrode in cylindrical 
geometry by solving Poisson’s equation for a case with only ions, represented by 
2000 charged rods, and showed that the virtual electrode is stable for long periods of 
time. In a simulation with both electrons and ions, however, Bames [16] found that 
electrons trapped in the virtual anode tend to wash out the anode potential, so that 
the virtual anode gradually disappears. Only in cases with very high ion currents did 
the virtual anode persist. He concluded that the electrostatic-inertial confinement 
device, as originally conceived, was unlikely to confine a fusion plasma adequately. 
Porter and Kievans [17] investigated the stability of electrons flowing among ions 
within a virtual anode, and concluded that the electrons are marginally stable to the 
two-stream instability. 

1.1.2. Multiple well hypothesis. Hirsch [18] observed high DT neutron yields (up to 
2 X 10 10 s -1 ) from a steady-state device using ion injection from six ion guns. These 
neutron yields exceeded the predictions of simple theories. One hypothesis to 
account for the observed neutron yields is the formation of multiple, concentric 
spherical potential wells of alternate sign inside the central sphere, with much higher 
circulating currents in the inside layers, as illustrated in figure 2. According to such a 
model, ions trapped in interior layers could have high energies, circulating currents 



Figure 2. Multiple virtual electrodes inside the accelerating grid. 
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Figure 3. The results of a neutron collimation study at -90 kV, 20 mA and p = IPa 
(deuterium) [18]. 


and densities, and could produce significant fusion energy yields. Hu and Kievans 
[19] modelled multiple wells for cases with zero angular momentum (perfectly radial 
motion). They calculated the attainable current amplification ratio (the ion current 
inside virtual anode/ion beam current outside virtual anode), and found that values 
as high as 150 were possible if large trapped ion populations could be confined in the 
outer region. 

The following observations appear to be consistent with the multiple potential 
well hypothesis: 

(i) high neutron yield; 

(ii) peaks in the curve of neutron yield versus radius measured by Hirsch [18] as 
shown in figure 3; 

(iii) peaks in the curve of electron bremsstrahlung yield versus radius measured 
by Hirsch [18] as shown in figure 4; 



Position with respect to cathode vertical axis, mm 

Figure 4. Electron bremssstrahlung yield versus radial position, as measured by Hirsch 
[18]. 
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Figure 5. Theoretical variation of dimensionless electrostatic potential Y versus dimen- 
sionless radius X, for a case with applied voltage = 500 V, current = 125 mA, grid 
gransparency = 90%, and very small spreads of energy and angular momentum. A 
double well is apparent. 


(iv) a large time-delayed puff of gas after beam injection cut-off measured by 
Hirsch [18]; 

(v) theoretical studies by Cherrington et al [20] for cases with injected beams 
having narrow spreads of energy and angular momentum show a slight double well, 
as illustrated in figure 5; 

(vi) electron beam probe studies by Swanson et al [21,22] of a spherical 
electrostatic confinement device at low pressures; 

(vii) peaks in the curve of electron density versus radius (figure 6) measured by 
laser heterodyne studies in a cylindrical electrostatic confinement device by Meeker 
and co-workers [23, 24]. This figure is for the case of electron injection; similar, less 
pronounced peaks were observed during ion injection. 


Electron Density 



Radius, mm 


Figure 6, Electron density versus radius measured by Meeker et al using a laser 
heterodyne interferometer during electron injection into a cylindrical grid with radius 
40 mm, at 16 A circulating current, p = 1 Pa (deuterium) [24]. 
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Other studies, however, have yielded results less indicative of multiple potential 
wells. Imel [25] and Black and Kievans [26] developed theoretical models of the 
electrostatic potential profile and electron and ion density profiles, using distribution 
functions tailored to various experimental conditions. They analysed three cases [27] 

(i) Ion beam injection at very low pressures (the Hirsch experiment). They 
concluded that multiple deep potential wells would probably not form under the 
conditions of the Hirsch experiment. 

(ii) Ion injection at higher pressures, where charge exchange produces a broad 
spread of ion energies (several experiments). They predicted a central plasma 
density close to the value measured in the Penn State experiment (2X 10 1S m -3 ), and 
found that substantial potential wells could be formed for p < 0.13 Pa (1 mTorr). 

(iii) Electron injection at low pressures (the Swanson et al and Black and 
Robinson experiments). They predicted deep potential wells for high perveance 
cases, as observed by Swanson ef al, and shallow potential wells for low perveance 
cases, as observed by Black and Robinson [28]. 

They did not report any cases with finite angular momentum and multiple potential 
wells. 

Gardner and colleagues [29] continued the ion beam experiment of Hirsch and 
used a microwave cavity to measure the number of free electrons in the plasma, but 
they were unable to reproduce the neutron yields observed by Hirsch. Baxter and 
Stewart [30] predicted nuetroa yields close to those observed in the Hirsch 
experiment by analysing ionization and charge exchange in the ion beam, without 
invoking the multiple potential well hypothesis. 

Peterson and Oleson [31] studied hemispherical focusing of a low-energy ion 
beam in a double-beam plasma device. When the ion beam energy exceeded T c the 
beam ions were slowed down by a potential hill (virtual anode) at the focus. They 
observed formation of a dense group of cold electrons (0.8 eV), which appeared to 
ride along with the ion beam, in addition to the background plasma electrons 
(4.5 eV) and primary electrons from filaments. Such cold electrons could tend to 
wash out potential wells. 

Nadler et al [32] built a collimated solid state proton detector to observe protons 
resulting from DD fusion reactions in a spherical electrostatic confinement device. 
The energy shifts produced by aluminium and lead foils verified that they were 
detecting 3MeV protons. At 12 mA cathode current, 30 kV voltage, in 0.5 Pa of 
deuterium fusion reactions between the focused ion beam and the background 
plasma were dominant. In future experiments at higher currents, the higher proton 
count rates should permit improved angular resolution, which could facilitate the 
detection of multiple potential wells. 

1.1.3. Regimes of operation. Miley et al [33] distinguished three regimes of ion 
injection into an inertial electrostatic confinement (IEC) device: 

(1) A discharge mode. In this mode a self-sustaining plasma discharge (like a 
glow discharge) generates ions throughout the chamber volume. Free electrons are 
produced by ionization and by secondary emission from grids. This mode requires 
comparatively high gas pressures (>1 Pa), which result in substantial charge 
exchange rates. 

(2) An emitter mode, based on a patent by Hirsch and Meeks [34]. In this mode 
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primary electrons from thermionic cathodes orbit in and out of an intermediate grid, 
ionizing background gas, at pressures typically 0.1-1 Pa. Ions produced by ionization 
are then accelerated into the central cathode. This mode was used in experiments by 
Meeker and co-workers [24, 35] to measure density profiles. Using this mode with 
deuterium gas, Miley et al [33] observed steady-state neutron yields of 10 s neutrons 
per second. 

(3) A beam injection mode. High-energy ion beams are focused into the central 
region, to produce a high-density virtual anode. This mode, which can operate at the 
lowest pressures, was used by Hirsch to obtain high neutron yields. This mode is 
most suitable for fusion reactors, because the low pressure minimizes ion energy loss 
by charge exchange. 

1.1.4. Reactor prospects. For steady state operation with radiative cooling of a 
tungsten grid intercepting 1% of the ion current and focusing the ions to a virtual 
anode at radius r 0 = O.Oli?, Lavrent’ev [36] found that a grid radius R > 10 m would 
be needed in order to achieve a net energy output. For the case of a pulsed reactor 
the required grid radius is given by 

R >0.Q01yt p <f> 3 A (r o /R)/4d (3) 

where y is the fraction of current intercepted by the grid, r p is the pulse length, d is 
the grid wire diameter and <j> A is the applied voltage (kV). For a case with y = 0.01, 
rjR = 0.01, <f> A = 100 kV, t p = 0.1 s and d = 0.001 m, this equation gives R >2.5 m. 

It is also possible to focus the injected ions temporally by increasing the voltage 
during injection [37]. Ions starting out later are accelerated more and catch up with 
the early ions, so that the ion density at the focus is much higher than the 
steady-state value for continuous injection. For example, starting with background 
plasma at radii between R = 1 m and 2 m with density n 0 = 10 1S m“ 3 , the focused 
density would be n = 6 X 10 25 m -3 . The converging ion beams could also be used to 
compress and ignite inertial confinement fusion target pellets. 

Although purely electrostatic plasma confinement might be attractive for a fusion 
reactor with a large radius or short pulse length, overheating of the grid wires is a 
serious concern. One way to reduce the heat load from charged particle bombard- 
ment of the grid wires is to shield them magnetically by passing high currents 
through them, as illustrated in figure 7. If electrons are injected between the grid 
wires, very few will contact the grid wires directly. Some electrons will gradually 
become trapped in the magnetic field and diffuse to the grid wires, but this process is 
orders of magnitude slower than direct bombardment. This magnetic shielding of 
grid wires, here called MEPC, has also evolved from open magnetic confinement 
systems. 



figure 7. Cutaway view of a few grid wires with currents flowing out (+) and into (— ) 
the plane of the drawing. The ellipses represent magnetic field lines. 
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1.2. Open magnetic confinement systems 

Minimum-B magnetic confinement systems offer good MHD stability, but suffer 
from rapid plasma loss along magnetic field lines [38]. Confinement times in simple 
magnetic mirror cells are typically of the order of ion-ion scattering times, and 
confinement times in magnetic cusps are tens of transit times, which typically makes 

Q s (fusion power)/ (input power) 1 (4) 

for ordinary mirrors and cusps. In order to attain Q > 10 without requiring very 
large size or extremely high magnetic induction, it is desirable to find a means for 
plugging end losses from magnetic mirrors and cusps. Many means of Q- 
enhancement for open magnetic confinement systems have been studied, including 
tandem mirrors, multiple mirrors, field reversed mirrors, Tormac, Surmac, plasma 
rotation, gasdynamic traps, radiofrequency plugging and electrostatic plugging 
[39,40]. Axisymmetric cusps have been suggested as electrostatic end plugs for 
tandem mirrors [41-43]. 

Several cusped magnetic fields have been considered for plasma confinement 
[44], as illustrated in figure 8. The toroidal multipole shown in figure 8 has N = 6 



Figure 8 . Magnetic field lines (smooth curves) in some cusped geometries. The symbols 
x and * represent current into and out of the drawing, respectively. 
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cusps (hexapole). In addition to these configurations, arrays of point cusps on the 
surfaces of polyhedrons have been studied by Keller and Jones [45] and by Sadowski 
[46-50]. All of these cusp geometries have been considered for electrostatically 
plugged devices. 

For plasma confinement by magnetic cusp fields, the loss area through a ring 
cusp is approximately 2jiR{28), where R is the ring cusp radius and 8 the effective 
half-width of the cusp gap through which plasma is lost. In a spindle cusp, the loss 
area through the point cusps is comparable, so the total loss area A loss = 8?r RS. The 
ion flux into the loss regions is approximately 0.25 nu, where u is the ion speed, so 
the confinement time in an unplugged magnetic spindle cusp device is approximately 

T«=nK/(0-25n«Ai oss ) = V /ItiuRS (5) 

where V is the plasma volume. This time is typically tens of bounce times back and 
forth inside the central field-free region. Spindle cusps have very small plasma 
volumes. If 8 ~p b then unfeasibly large R and B would be required for a reactor. If 
8 ~ (p„Pi) 1/2 (a hybrid gyroradius), as indicated in some cusp experiments and in a 
computer simulation [51-53], then a reactor using a ‘picket fence* cusp (such as the 
spherical cusp and linear set of ring cusps in figure 8) is feasible, but high magnetic 
fields and large radii are needed [54]. 

For a linear set of ring cusps (figure 8) with distance Lx between the ring cusps, 
the plasma volume between cusps is V = lar^Ly, where r p is the effective plasma 
radius, and the loss area A, oss ~ 4? zR8. For this case 

t — rpLJuRS. ( 6 ) 

This equation, which estimates the plasma loss time by free-flow-out ring cusps with 
an effective loss width 25 in each ring cusp, will be used later to estimate the 
effectiveness of electrostatic plugging. 

Lavrent’ev and colleagues [5,55] proposed electrostatic plugging of cusps to 
reduce the loss rate of plasma flowing along magnetic field lines out of the cusp gaps, 
and conducted early experimental studies of this technique. If the plasma potential is 
highly negative relative to the walls, the ions will be confined electrostatically. Only 
electrons will be energetically able to pass through the narrow cusp gaps (plus a few 
ions in the high-energy tail of the Maxwellian), so the width of the untrapped 
electron stream flowing through the cusp gap may be reduced to 8 ~ p e . Cathodes 
outside the cusp gap reflect escaping electrons back into the plasma. Such 
electrostatically plugged cusps are equivalent to magnetic shielding of the grid wires 
of an electrostatic-inertial confinement device. In electrostatically plugged cusps, the 
ions no longer flow freely out of orifices, so the ffee-flow confinement time estimates 
do not apply. Instead, the ions are confined electrostatically by the potential well of 
electron space charge, and electron confinement is limited by diffusion across the 
magnetic field. 

Moir discussed the use of electrodes to contain warm plasma for stabilizing open 
magnetic confinement systems, such as a toroidal quadrupole cusp [56]. The 
augmentation of toroidal magnetic confinement by strong electrostatic fields was 
studied by Daugherty et al [57] and by Stix [58-60]. Jones [61, 62] suggested the use 
of electron or ion injection to create radial electrostatic potential variations, instead 
of using plugging electrodes. 
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Figure 9. An electrostatically plugged toroidal quadrupole magnetic cusp and the 
resultant distribution of electrostatic potential <j>. The broken line represents the 
potential with no plasma present. 


2. Theory 


2.1. Electrostatic potential 

Magnetic electrostatic plasma confinement, which may be accomplished with any 
magnetic cusp configuration, is illustrated in figure 9 for the case of a toroidal 
quadrupole cusp. High-voltage electrode rings are placed in each of the cusp gaps, as 
shown. The anodes are biased positive, and the cathodes negative. Operation of the 
confinement system depends mainly upon the voltage <£ A applied between the 
cathode and anode, and little upon the relative location of the ‘ground’ potential. 
From figure 9 it is evident that 


4>A - <£e + 


( 7 ) 


This potential shape is the same as that of a simple tandem mirror. 

In a vacuum, the interior of the device will be near the anode potential (broken 
curve of figure 9). Plasma may be produced by electron beam injection into 
low-pressure gas, by rf (radiofrequency) heating, by plasma gun injection, by 
laser-pellet heating, by neutral beam injection or by other means. Electrons exiting 
the cusps are reflected by the negative voltage of the cathode, but the exiting ions 
are not initially confined. Since electrons are lost less rapidly than ions, the plasma 
develops a slight negative charge, and the plasma potential becomes negative 
relative to the anodes (the full curve of figure 9). If the anode gaps are narrow 
enough that their potential is not entirely shielded out by the plasma in the anode 
regions, then the potential there will be near the anode potential and higher than 
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Figure 10* Hypothetical cqui potential surfaces in the anode region, showing saddle- 
shaped electrostatic potential distribution. 


the plasma potential, forming a potential hill fa for the remaining ions. Only a slight 
fractional charge imbalance is required to set up a potential hill many kV high. Then 
only ions with kinetic energies greater than 6(f), can escape along the magnetic field 
lines out of the cusps. The two-dimensional potential distribution (figure 10) is 
saddle-shaped in the anode regions, with the saddle point lying an amount A<£ below 
the anode potential. 

The potential hill for ions approaching the walls (at anode potential) is even 
higher, so most ions will go out through the cusp gaps as soon as they acquire 
enough energy to overcome the barrier fa; no ions acquire enough energy 
collisionally to surmount the barrier (<£j + A<£) and reach the walls. The ions are 
electrostatically confined in a negative electrostatic potential well produced by a 
slight charge imbalance. The ion Larmor radius is less significant than usual, because 
the ions are electrostatically reflected in the boundary layer. The maintenance of a 
potential difference along the magnetic field lines is facilitated by the boundary 
condition that the walls are very close to the plasma only in the anode gaps. Some 
energetic ions are replaced by cold ions as a result of charge exchange with neutral 
gas. 

The central plasma region is uniform in density and temperature, free of electric 
and magnetic fields, and surrounded by a thin edge layer. The outside boundary of 
the edge layer is the magnetic surface along which electron orbits graze the anodes, 
and the inside boundary is the uniform central plasma region. The boundary layer 
thickness at any location with magnetic induction B is found by magnetic flux 
conservation to be 

A r**a(RBJrB), (8) 

where a is the anode gap half-width, R is the ring cusp anode radius and r the radius 
of the boundary layer. If B a - 5T, B = IT, R = 4 m, r = 3m and a ■ 1.5 mm then 
Ar = 10 mm. This thickness is much less than the plasma size, so the plasma has a 
‘sharp boundary’. 



Magnetic electrostatic plasma confinement 


1553 



Figure 11, Electron orbits near a magnetic cusp. Smooth arcs represent magnetic field 
lines. Dots represent ions. Electron trajectories are as follows: 1 = reflected off plasma 
boundary, 2 = passing through anodes and reflected by cathode, 3 = magnetically 
trapped in boundary layer, 4 = electrostatically trapped by anode potential well, 

2.2. Electron orbits 

There are several types of electron orbits, illustrated in figure 11: 

(1) Electrons that stay in the central plasma and are reflected geometrically from 
the convex edge layer. 

(2) Electrons from the central plasma that can pass through the anode gap 
regions and be reflected by the cathodes. 

(3) Electrons that are magnetically trapped in the edge layer. (Many of these are 
mirror-reflected away from the cusp gaps.) 

(4) Cooler electrons that are electrostatically trapped by the anode potential. 
Most of these electrons are produced by ionization of incident neutral gas atoms 
along the slope of the potential well. These trapped electrons are deleterious, 
because they increase the self-shielding voltage drop 

Electrons may be produced by thermionic emission from an electron gun, by 
secondary emission from the plugging cathodes and by ionization of neutral gas. 
Much of the injected electron beam is reflected by the plasma and lost to the 
cathodes. This portion of the injected beam does not contribute to plasma density 
build-up. Yushmanov [63] found that only a small fraction of the injected beam 
current is effectively trapped. The time required for plasma density build-up by 
electron beam injection (many particle confinement times) exceeds the duration of 
the flat-top magnetic field in some pulsed experiments. 

Many of the electrons produced by ionization will be magentically trapped in the 
edge region. The untrapped electrons will circulate back and forth through the 
plasma, occasionally passing out through the anodes and being reflected by the 
potential barrier </> e . Coulomb collisions and electric field fluctuations cause diffusion 
in velocity space. Trapped electrons can become heated and detrapped, and free 
electrons can be deflected and trapped either magnetically or electrostatically [64]. 

There are two main electron loss processes: trapping by the magnetic field with 
subsequent diffusion to the anodes, and collisional diffusion in velocity space over 
the potential barrier <j> e with loss to the cathodes. (Loss to the cathodes converts 
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electron kinetic energy into electricity, returning energy to the power supply.) 
Recombination with ions is negligible. 

2.5. Loss rates over potential barriers 

The loss rates along the magnetic field by velocity-space diffusion over the potential 
barriers have been calculated by Cohen et al [65], based on the method of Pastukhov 
[66]. For a plasma with singly charged ions, the resulting particle loss rates along the 
magnetic field are 

(dn k /dt\ = -4 n k Z k exp(-y k )I(l/y k )/[jc m y k z k G(RZ k )] (9) 

G(RZ k ) = (1 + 1 IRZ k ) m ln{[(l + MRZ k ) m + 1]/[(1 + 1 /RZ k )) m - 1]} (10) 

10-bk) = 1 + 0-5(n/y k ) m exp(y*)[l - erf(y^)] 

= 1 + 1/(2 y k ) - 1/ (4y*) + 3/(8y 3 *) + . . . (11) 

where the subscript k denotes either electrons or ions, n k is the particle density, 
Z k = 1 for the electrons and 1/2 for ions, R is the magnetic mirror ratio, v k is the 
characteristic collision time for electrons or ions, and y k = z k e<f> k lT k where z k is the 
charge number (= 1 for singly charged ions, = -1 for electrons), e is the electronic 
charge and <f> k is the potential barrier for eectrons or ions (here z e and </> e are taken 
to be positive quantities). 

The energy loss rates along the field are given by 

1.5[d(n A T,)/dt] l( = z k e4> k [\II(lly k ) + 1.5/y*](dn*/d0„ (12) 

A similar estimate of loss rates was made by Sizonenko and Stepanov [67]. These 
equations give the rate at which particles or energy are lost by electron or ion 
diffusion in velocity space over their potential barriers </> e and In view of the 
exp(-_y ; .) factor in the particle loss equation, the confinement time should increase 
roughly expontially with y k . Computer simulations by McHarg and Oakes [68, 69] 
show such an exponential increase of parallel confinement time with <t> k /T k . In the 
limit of high plasma density or low temperature, collisions may be so frequent that 
these ‘Pastukhov’ equations lose validity, and the equations of Rognlien and Culter 
[70] should be used. In MEPC devices the dominant electron losses are by diffusion 
across the magnetic field, but ion confinement is almost completely electrostatic. 

Parks and Sleeper [71] derived an equation predicting ion end loss rates in the 
transition region between the non-adiabatic region (central plasma) and the 
adiabatic region (boundary layer). In the non-adiabatic limit, their equation agrees 
with the result of Yushmanov [72] and in the adiabatic limit, their equation is similar 
to the results of Chemin and Rosenbluth [73], but with loss rates a factor of 3-4 
higher, due to effects of magnetic field curvature. 


2.4. Self-shielding 

The uncompensated electronic space charge tends to shield out the anode potential 
by an amount A </>. Assuming azimuthal symmetry, the electrostatic potential 
distribution can be calculated from a two-dimensional Poisson equation in the anode 

gap 


(l/r)(d/dr)(rd4>/dr) + d 2 <j>ldz 2 = ( n - n{)e/e 0 . 


(13) 
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Due to the potential hill c& ; . only a few ions are passing through the anode gaps, so 
n ia «n, as will be discussed later. With comb probe measurements in some 
experiments it was found that the distribution of electrons in the anode regions had 
bell-shaped distributions across the magnetic field. In order to model a variety of 
conditions, a Lorentzian-shaped density profile 

n=nj(l + z 2 /w z ) (14) 

is assumed, where n a is the electron density at centre of the gap and the parameter w 
characterizes the half-width of the distribution. Assuming this electron density 
distribution in the anode region, with « a = 1.5 X 10 17 m~ 3 , a - 2 mm and w = 0.2 mm, 
results in the electrostatic potential distribution of figure 10. The depth of the saddle 
point relative to the anodes is called A</>. If the anode length along the magnetic field 
h»a, then a one-dimensional equation 

d 2 <j>ldz 2 = ( n - n t )e/ s 0 (15) 

can provide an accurate calculation of the potential distribution <f>(z) across the 
magnetic field. Assuming 4>(-a) = 4>(a) = 0, solution of this equation yields 

A(j> = (jm^ewa f2e 0 )[(2/ jr) tan ~*(a/w) - (wjna) In(l + a 2 [w 2 )\. (16) 

Since high-energy electrons flying through the anode gap have trajectories extending 
one or two Larmor radii on each side of the centre, it is assumed that w ~ 2p„. The 
resulting values of A 4> are shown in figure 12 as a function of n & and B , assuming 
7^=20keV and a = 2 mm. If n a = 10 19 m -3 and B = 6T, then Ac£ = 74kV. These 
could be typical parameters for an MEPC reactor. A comparable result was found 
by Ware and Faulkner [74], assuming a triangular density distribution. The quantity 
in brackets varies from 0.96 at w/a = 0.01 to 0.67 at wfa = 0.2, so the approximation 

A(f> » 7m a ewa/2£ 0 (17) 

gives a slight overestimate of A<j>. Expressing T ek in keV and A<p in kV, this becomes 

A <£ = 6 x 10" 15 T lfn a a/B (18) 

which yields A<f> ~ 89 kV (20% high) for the above example case. The maximum 
electric field occurs at the anode wall 

£max = -(d<A/dz) ;=fl = (n a ew/£o) arctan(a/w) ~ A <fla. (19) 

For typical values of MEPC experiments, the electric field energy density 
(l/2)£ 0 E 2 m ^«B 2 /2» 0 . 



Figure 12. Variation of self-shielding potential sag with magnetic induction, for various 
anode gap peak electron densities n a1 assuming w = 2 p c , a = 2 mm, T = 20 ke V. 
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2.5. Electron density in aode gaps 

Five phenomena affect the ratio njn: ions in the gap, magnetic reflection, 
electrostatic acceleration, cold trapped electrons and diocotron oscillations. 

2.5.1. Ions in the gap. Assuming a quasicylindrical plasma with radius r v and length 
L, the rate at which ions are lost is 

7 tr 2 pLn /t ; = N27tR28n ia u iSL (20) 

where n is the central plasma density, is the ion confinement time, N is the 
number of ring cusps, R is the ring cusp radius, S is the effective half-width of cusp 
gaps for ions, n; a is the density of ions flowing over the barrier and u ia is the speed at 
which ions pass over the barrier. Ions just barely leaking over the potential barrier 
will have low velocities « ia ~3xl0 4 ms _1 , and they will be restricted to a narrow 
corridor of half-width S~w. At equilibrium the ion particle confinement time will 
equal the electron confinement time, which is roughly equal to the diffusion time t diff 
across the magnetic field. Thus, the ion density in the anode gaps is estimated to be 

njjn = r^Li/4Rwu ia r diS (21) 

where L x = distance between ring cusps. For the example MEPC reactor parameters 
of table 1, assuming L x = 1.5 m, it is estimated that n jn ~ 0.02. 

2.5.2. Magnetic reflection. Sidorkin and Lavrent’ev [75] studied electron reflection 
by multipole cusp gaps, as illustrated in figure 13. They showed that electrons with 
a > a crit are reflected by the cusp field, where 

sin(a crit ) = [(1 + + x 2m )x‘S^] m (22) 

*o = [(m - 1 )/(m + l)] 1/2m (23) 

x = r/R, r is the initial electron radius, R is the cusp gap radius and m is the number 
of pairs of multipole conductors = N/2. The values of a crit are plotted as functions 


Table h Example MEPC reactor parameters* The plasma parameters were estimated 
using approximate scalings of attainable density, temperature, confinement time (section 
3.7.1) and Q (section 4.1). 


Ring cusp magnetic induction, B a — 6 T 

Applied voltage, <j> A — 400 kV 

Plasma radius, r p = 3 m 

Ring cusp anode radius, R = 4 m 

Ring cusp anode gap half-width, a =2 mm 

Estimated, parameters 
Central plasma density, n ** 10 20 m -3 
Electron and ion temperatures, T c , 7^ « 20 keV 
Magnetic field at plasma boundary, 1.3 T 
Electron Larmor radius in anode gap, ** 8 x 10 -5 m 
Electron Larmor radius at the pLasma boundary, 
p b —3.8 X 10 -4 m 

Self-shielding voltage drop, A <j > » 80 kW 
Anode gap peak electron density, 1.3 X 10 19 
Energy confinement time, t E± » 6 s 
Power gain ratio, Q « 10 


Equations 

(57) and Jupiter-2M 
(48), (49), (58) 
after (37) 
after (37) 
after (37) 


(47) 

(16) 

(40) 

(60) 
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Figure 13. Current-carrying conductors of a multipole cusp (left) and definition of the 
electron trajectory angle of incidence a relative to the midplane of the cusp gap. In this 
figure N = 8, m = 4. 

of initial radii in figure 14, for various values of m. In multipoles with large m the 
electrons have smaller critical angles, so more of them are magnetically reflected. 

Electrons starting off at r/R = 0.5 near an octupole cusp (m = 4) have x 0 = 0.938 
and a crit = 29°. Experimentally measured electron beam transmission is consistent 
with such theoretical predictions. Assuming that the central plasma electrons have 
isotropic distributions, the total loss fraction through N cusps from a given location 
is approximately 

N 

total loss fraction 0.5 2 [1 ~ cos(a crit )/]. (24) 

/=i 

If rt p = density of electrons that are able to penetrate into the magnetic cusps, the 
ratio n p /n can be estimated from the volume-averaged value of this total loss 
fraction, but this calculation has not yet been performed. 

2.5.3. Electrostatic acceleration. Electrons that succeed in entering the cusp gaps are 
accelerated by the positive applied voltage there. Moir et al [76] showed that the 



Figure 14, Variation of critical reflection angle with radial position from which the 
electron trajectory starts, for various multipole orders. Smooth curves are theoretical, 
and full points are experimental data for m - 4. 
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electron density in the gaps is lower than the magnetically penetrating density n p by 
a factor 

njn p = exp(y) erfc(y 1 ' 2 ) = (ary) -1/2 (l - l/2y + 3/4 y 2 - 1 5/8y 3 + . . .) (25) 

where y = <f>JT e , and erfc is the complementary error function [77]. For values of 
(f>JT e . = 3 to 6, this ratio is rtjn p = 0.29 to 0.22. 

2.5.4. Cold trapped electrons. Magnetic reflection and electrostatic acceleration tend 
to reduce njn, but electrons produced by ionization of neutral gas tend to become 
electrostatically trapped in the anode region (orbit 4, figure 11), increasing njn. If 
classical cross-field diffusion were the only loss process, then the accumulation of 
trapped electrons could increase n a and A <f> to unacceptably large values, but the 
diocotron instability may help remove cold trapped electrons without seriously 
impairing hot electron confinement. 

2.5.5. Diocotron oscillations. The electric field due to electron space charge is zero 
at the middle of a ring cusp gap, increasing on each side. The resultant E X B drift 
velocity varies spatially, and the shear of this velocity gives rise to the diocotron 
(slip-stream) instability [78,79]. Computer simulations by Levy and Hockney [80] 
showed good agreement with theoretical instability growth rates. Pankrat’ev et al 

[81] found that the long-wavelength instability (A > lOw, where w is the electron 
stream half-width) occurs when the electron drift velocity is near the wave phase 
velocity, and it is suppressed by a nearby conducting wall. However, Gordienko et al 

[82] found that the presence of two groups of electrons (trapped and passing) in the 
anode gaps broadens the unstable range of wavelengths and makes diocotron 
oscillations unstable whenever there is a vacuum gap between the electron stream 
and the anode wall. 

The short-wavelength instability (A < 0.2 a) occurs near the electron cyclotron 
frequency, and has a growth rate 

r/ &*ce ~ (<? /2) exp(-2lq) (26) 

where q = <o p J coZa- These diocotron oscillations play a beneficial role by removing 
cold trapped electrons from the anode region, where they would broaden the width 
w of the electron density distribution and increase A </>. The instability growth rate is 
small if q < 0.2, which is equivalent to the condition that 

rt a *£ (2 X 10 18 m -3 )i?| (27) 

where B s is the magnetic induction in the gap (T). Thus, as n a decreases well below 
this limiting value, the short-wavelength diocotron instability turns off. At B a = 5 T, 
this limit is n a ^ 5 X 10 19 m -3 , which is less restrictive than the value of « a =s 
1.4 X 10 19 m -3 required to keep A <j> < 100 kV (figure 12). 

Yushmanov [83] studied the build-up and removal of trapped electrons in the 
anode regions, and determined the required anode dimensions. The width h of the 
anode (figure 11) should be in the range 

5 a^h^ h max m 6.7 (aelN){n t /n a )(p 2 Ja z )(BJB b )R (28) 

where a ~ 4.7, e = trapped electron pumpout rate/diffusion rate = 1/3, N is the 
number of cusp gaps and njn a is die trapped electron density/anode electron 
density = 1/3. Taking N = 6 in the example case of table 1, it is found that 
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Figure IS. Maxwellian distribution of central plasma electrons (top): shifted Maxwellian 
distribution in anode region, assuming ^ =2 T e (middle); relative diffusion rate in anode 
region, with effect of diocotron oscillations on cold electrons (bottom). 


l^niax ~ 4 cm. Low-frequency oscillations can selectively pump out the trapped 
electrons without seriously increasing the loss rates of the free plasma electrons. 

If the electrons in the central plasma have a quasi-Maxwellian distribution, the 
distribution will be upshifted in energy by an amount ed>, as the electrons enter the 
anode regions, as illustrated in figure 15. Electrons produced by ionization in the 
plasma edge region will have a spread of energies. The diffusion rate is high for 
low-energy electrons, due to interaction with diocotron oscillations, and it may be 
nearly classical ( D « W~ m ) for high-energy electrons. The resulting electron energy 
distribution in the anode regions is expected to resemble the shifted Maxwellian 
distribution of figure 15 (middle). 

Some magnetically trapped electrons have low enough energies that they can 
immediately interact with diocotron oscillations; but most of them, having been 
accelerated by the anodes, must be decelerated by collisions before the diocotron 
oscillations interact strongly with them. For one example case [84] the trapped 
electron fraction was estimated to be nJn. A ~ 0.24. Considering ions, magnetic 
reflection, acceleration, cold trapped electrons and diocotron oscillations, there will 
probably be about an order of magnitude density ratio between the anode region 
and the central plasma 

njn — 0.1. (29) 

The ratio would probably be lowest for cases with a large volume of field-free 
plasma, narrow anode gaps and low neutral gas pressure; and it could be hear unity 
for cases with spindle cusp magnetic field, wide anode gaps and high neutral gas 
pressure. Values inferred from data in a variety of experimental conditions range 
from 0.01 to 1. 
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2.6. Plasma potential 

At T e >l keV and «<10 2O m -3 recombination is negligible. Ions are confined 
electrostatically with energy loss by charge exchange and by diffusion in velocity 
space over the potential barrier. The ion and electron particle conservation 
equations may be written 

drtj/df = 5; - (dn/dOn (30) 

dn/dt = Si + I/eV - (dn/dt\-n/t L (31) 

where S K is the ionization source, I is the electron current from cathodes into the 
plasma and r ± is the electron loss time by diffusion across the magnetic field. At 
equilibrium the plasma potential will adjust itself so that these electron and ion loss 
rates are equal [85] 

IleV - (d n /dr), — n/r x + (dn E / dr),, = 0. (32) 

Since 

= <#> A - - A 4> (33) 

these two equations can be solved simultaneously [86] for the two unknowns <p e and 
<f > i5 using equation (9) for the parallel loss rates. The plasma potential distribution 
has been computed self-consistently in some two-dimensional models. 

2.7. Two-dimensional models 

Brunei et al [87] solved a two-dimensional (r,z) Poisson equation, assuming nearly 
Maxwellian distribution functions dependent on the constants of motion. They 
calculated the distribution of potential and particle densities in an electrostatically 
plugged spindle cusp corresponding to the KEMP device. They found that the 
electron and ion densities were both maximum on-axis. The electron density had a 
secondary peak in the ring cusp anode, and the ion density was very low in the 
anode regions. Matte and Lafrance [88] developed an improved numerical method 
for solution of the Poisson equation in MEPC devices. 

Shoucri et al [89] simulated an MEPC plasma with about 5000 electrons and ions 
(M/m = 16, where M is the ion mass and m the electron mass) in a 128x50 
cylindrical grid, using a capacitive matrix Poisson solver, with parameters similar to 
those of the KEMP II experiment. The simulation showed the maintenance of a 
two-dimensional electrostatic potential well, with saddle points in the anode gaps, as 
expected. Electrons were very well confined, but ion confinement was limited by 
plasma self-shielding in the point cusp gaps, especially during electron injection. The 
confinement time observed in the simulation was scaled up by the square root of the 
experiment/simulation density ratio to yield a predicted confinement time of 
50-75 ps for the KEMP II experiment, which was consistent with experimental 
measurements. 

2.8. Energy confinement time 

The non-radiative cross-field electron energy loss time may be written 

1/^EJ. — 1/tcondri" 1/ + ^trap) (34) 

where T cond is the time for heat loss by conduction ~T.5T diff , z dia is the characteristic 
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time for diffusion across the magnetic field and r trap is the time for free electrons in 
the central plasma to become trapped by the boundary magnetic field. 

One estimate of the trapping time is [90] 

Ttrap = 8ne 0 r p (2mvB b /ne 3 In A) 1/2 (35) 

where v is the electron velocity and In A is the Coulomb logarithm. For reactor 
conditions T trap « T di(f , and 

Tej. “O-tfrdiff. (36) 

Nearly classical transport has been observed in some MEPC experiments, while 
rapid anomalous transport has been seen in others. 

2.8.1. Classical transport. Pastukhov [91] estimated the classical diffusion time to be 

^diff — 2r C j<?K / Sap b p 3 (37) 

where z ei is the electron-ion momentum-transfer collision time, V is the plasma 
volume, S is the plasma surface area, p b = (2 mkT e ) l/z (eB b is the electron Larmor 
radius at plasma boundary with B b = [2 p 0 nk(T e + 7])] 1/2 , p a = (2mkT e ) m leB a is the 
electron Larmor radius in the anode gap and a = a parameter varying from 1.33 at 
Ti/T e = 0 to 4.6 at 7]/7^ = l. A similar estimate was derived by Dolan et al [92]. 
Assuming In A = 16 in a hydrogenic plasma, the electron-ion momentum-transfer 
collision time [44, 93] may be reduced to the form 

T e j = 9.4 X 10 u T 3 £/n (38) 

where T ek is in keV. For spindle cusps, V/SR -0.01-0.03, depending upon 
dimensions, plasma heating and confinement effectiveness. (Low-pressure plasma 
have smaller volumes.) This low ratio plus diffusion enhancement by electron 
collisions with neutral gas give spindle cusps poor confinement times. 

Assuming classical diffusion in a spindle cusp with radius R and length L, 
Lavrent’ev [94] derived analytical expressions for the electron density profile, the 
plasma volume, the average electron density and the electron confinement time. The 
electron energy loss rate is 

P e = [4>, + A4> + 2T e (cosh(p) - 1 )]7tR 4 D 0 n 0 p/rlL sinh(p) (39) 

where D 0 is the collisional (classical) diffusion coefficient, n 0 is the central plasma 
density, p = (Tj./T],) 1/2 and r a is the point cusp anode radius. This power loss is 
minimized when p ~ 5. Lavrent’ev found that the edge electric field can cause 
significant increases of electron transport rates in. small devices with substantial 
neutral gas pressures. For the parameters of the Jupiter- 1M device the electric field 
and electron mobility should result in about a 45% increase in the electron 
cross-field transport above the classical diffusion rate. In the case of a reactor with 
‘burn-out’ of the neutral gas, however, the electron mobility would be small, and the 
electric field effect would be diminished. 

Lavrent’ev [95] showed that angular momentum from EXB rotation produces 
an effective potential well <p s that greatly inhibits axial ion losses, even if the point 
cusp electrostatic barrier 4>i = 0. A similar effect acts on electrons, but is injM times 
weaker. Near the ring cusp the electric field effects dominate the ion motion. 
Dolgopolov et al [96] computed ion trajectories in an electrostatically plugged 
spindle cusp, and found that the magnetic field gradient force and the centrifugal 
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force (related to the azimuthal EXB drift) tend to keep ions from escaping through 
the point cusps. This inhibition of point cusp ion losses makes plasma confinement 
effective in a linear set of ring cusps. 

For a linear set of ring cusps , the volume-to-surface area ratio V / SR =* 0. 1-0.3 is 
an order of magnitude better than or spindle cusps. Assuming V/ SR =0.02 and 
T e = Ti, equation (37) becomes 

T diff * 0.09r ei ai?/p b p a . (40) 

For the example MEPC reactor parameters case (table 1) it is found that 
T ei = 0.84 ms, the diffusion time r difi — 20 s and t Ejl ~ 12,s. Since experiments have 
attained about half the classical value, a value of t ej . ~ 6 s is listed in table 1. This 
near-classical transport would make an attractive reactor, but anomalous transport 
could prevent a reactor from being successful. 

2.8.1. Anomalous transport. Pastukhov and IPgisonis [97,98] studied anomalous 
electron transport in the boundary layer of an MEPC plasma. When 7] « T e , the ion 
acoustic instability is dominant, giving rise to rapid electron diffusion through the 
plasma edge layer. The long-wavelength ion-acoustic instability can occur at fcp e ~ 1 
and to ~ ((a ce (Od ) 112 when the electron density-gradient drift velocity v Vn >c s . The 
corresponding maximum diffusion coefficient jD max ~c s p e , which has Bohm-like 
scaling. The resulting anomalous transport tends to broaden the electron density 
profile to a condition where v Vtt ~ c s . The resulting theoretical electron density 
profile is given by 

n(r) = n 0 exp(-ec s A/cT e ) (41) 

where n 0 is the central plasma density and A is the magnitude of the magnetic vector 
potential. The waves are partially absorbed by the ions, and may heat them to the 
point where the instability turns off. When the ion acoustic instability 

becomes a Buneman instability, and the instability growth rate decreases by a factor 
of about (m/M) yz . Once electrons diffuse outside the bounding flux surface (the 
surface which just touches the anodes), they are magnetically mirror-confined (orbit 
3 of figure 11) until they scatter into the loss cone and travel along the magnetic field 
to the anodes. For a linear quadrupole with a/p e .<2(Mlm) in , the anomalous 
electron confinement time is 

t = {0.4 + 0.3 exp[2(a /p e )(m / M ) 12 ]}t 0 (42) 

where t 0 is the electron collision time. This case is similar to that of the ATOLL 
experiment, and the resulting confinement time is much shorter than with classical 
diffusion. Anomalous transport would spoil the energy balance in a reactor, resulting 
in a low fusion power density. 

2.9. Energy balance 

As neutral atoms impinge on the plasma some are ionized in the boundary layer. 
The resultant ions are accelerated as they fall down the potential hill into the plasma 
region. However, Pastukhov [91] found that this ion heating process is inefficient, 
and that ion heating by electrons is very slow. Ions are also heated by plasma waves. 
Ion energy is lost by charge exchange and by convection, as heated ions escape over 
the potential hill <f> : . Since ions with energies Wi > efa are lost out the cusps, their 
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distribution is truncated at higher energies, and the resulting fusion reaction rate is 
reduced [99]. 

Electron energy may be supplied by electron beam injection, neutral beam 
injection, wave heating and heating in the edge electric field. Electron energy is lost 
by electron flow over the potential barrier to the cathodes, by transfer to waves, by 
thermal conduction and convection across the magnetic field and by inelastic 
collisions with ions, atoms and molecules. 

Approximate energy balance equations for the central plasma electrons and ions 
may be written 

1.5 d(nT & )/dt = I4>JeV + P E + 0.2 5nf{<rv)W a f ae - 1.5[d(n7;)/dfJ l( - 1.5nTJz E± 

- Prad - 1 -5n(T* - 70/T eq - S m Wioss - Pw (43) 

1.5 d(«7])/dr = 1.5ft (7^, - 7])/-r eq +fiS m W> + p + 0.25nf{ov)(f ai W a - 1.8 efr) 

- 1.571(1 ~fi)S m - 1.5(dft,7j/dr)n (44) 

where P e is that portion of electron heating by the edge electric field which affects 
the central plasma, r eq is the electron-ion equipartition time, P rad is the radiative 
power loss, W loss is the average electron energy loss per incident neutral atom due to 
ionization, P w is the electron power loss to plasma waves, is the fraction of incident 
neutrals which are ionized (instead of causing charge exchange), f ae and f ai are the 
fractions of fusion product alpha energy W a which are transferred to electrons and 
ions, S m is the volume-averaged neutral input rate, VV, is the average energy acquired 
along the slope of the potential well by ions, ( ov ) is the DT fusion reaction rate 
parameter and p is the ion heating by plasma waves (such as ion acoustic). It is 
assumed that the average energy of each fuel ion undergoing a fusion reaction is 
about 0.9e<£j. Since the alpha particles are poorly confined, f at « 1 and f ai « 1. 
Additional terms would be needed for neutral beam injection, microwave heating, 
etc. These equations are nonlinear, because several parameters depend on the 
density and temperatures. Similar equations have been solved numerically for 
several cases [84]. Typical conditions satisfying the plasma potential and energy 
balance equations were found to be 


a ~ 0.5 
<£i/<£a~0.3 
A<£/<j> A «0.2 
Ti *= 0.05<£ a 
2>0.05<£ a . 


(45) 

(46) 

(47) 

(48) 

(49) 


In order to confine plasmas with fusion reactor parameters, applied voltages 
<£ a ^300 kV are desirable [100]. 

Karpukhin et al [101] studied particle and energy balance in an electromagnetic 
trap with r p ~0.3m, 5 a = 1-3 T, <t> A ~ 25-50 kV, V ~ 0.5 m 3 , a = 1.5 mm and an 
electron injection power of 500 kW. They predicted values of n ~ 3 X 10 19 m" 3 , 
t e ~ 0.1 -0.2 s and 7|~ 0.5-1 keV. For an octupole electromagnetic trap they found 
the optimum electrode gap spacing a — 0.5-1 mm, which maximizes the injection 
efficiency, plasma density and confinement time [102]. Energy balance calculations 
will be discussed further in section 4.1. In addition to spindle cusps, multipoles and 
ring cusps, two additional MEPC configurations have been studied recently: the 
Polywell tm concept and the modified Penning trap. 
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Figure 16. Electron injection through multiple point cusps in a Polywell 1 " 1 device. The 
full curves represent magnetic field lines, which are pushed outward by electron 
diamagnetism during higb-bela operation [109]. 


2 JO. The PolywelT m concept 

Bussard [103] proposed electrostatic confinement inside a polyhedral array of 
magnetic point cusps (figure 16), called a ‘Polywell Lm \ with electron injection at 
10-100 kV along point cusp axes to provide a quasi-spherical negative potential well 
capable of trapping and confining ions. The ions are to be injected at very low 
energies (<100 eV) from the system boundary (r = R), along cusp axes or the axes 
of the polyhedron vertices. The electron distribution in the ‘mantle’ region (r > r c ) 
outside the core is essentially spatially isotropic, but the electrons do not have a 
Maxwellian velocity distribution. The ion flow is expected to remain spatially 
anisotropic in the mantle region and to retain its non-equilibrium (nearly mono- 
energetic) nature over typical ion lifetimes. 

The degree of convergence of ions flowing radially to the central core is 
determined by angular momentum conservation to be 

rjR = (Wje^ m (50) 

where W { is the ion energy perpendicular to the radial direction and </>, is the 
accelerating potential [104], The PolyweH mi concept is based on high-beta electron 
confinement by the magnetic field. Watrous et al [105] showed that the required 
diamagnetic currents could be achieved with reasonable electron injection condi- 
tions, yielding a relatively field-free internal region that enhances electron confine- 
ment. Krall et al [106] also studied the density and pressure anisotropy profiles 
resulting from various ion sources, and assessed processes which may degrade ion 
focusing, such as angular deflection and energy upscattering. The concluded that 
isotropy of the plasma core inhibits the Weibel instability. 

Solving the Poisson equation for spherically convergent charged particle flows. 
King and Bussard [107] found that, with a high degree of electron recirculation 
within the system (with long electron lifetimes), central plasma densities of about 
10 2s m~ 3 could be attained at high electron gun currents. Bussard and King [108] 
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distinguished two bounding operating regimes: a low-electron-beta mirror regime, 
where magnetic trapping and reflection are dominant, and a ‘whiffle-ball (WB) 
regime', where a high-electron-beta plasma pushes the magnetic field outwards, so 
that the plasma is situated inside a quasi-spherical magnetic surface with small holes 
corresponding to the point cusps. In the WB regime they find that the electron 
recirculation factor <2 rO (representative of the central plasma electron lifetime) 
increases almost linearly with n c rl, where n c is the central electron and ion density 
within the core radius r z . They advocated plasma start-up at low magnetic field, in 
order to avoid the mirror regime, which requires higher electron beam currents. 

Three issues arise with the Polywell tm concept: 

(1) It may be difficult to achieve a high degree of ion focusing (small W x ), due to 
ion reflection by the scalloped electrostatic potential well boundary [109] and ion 
deflection by wave-particle interactions. Larger W ± would result in a larger core 
radius r c and reduced central plasma density n c . Further research is needed to clarify 
the attainable degree of focusing. Even if great focusing were not achieved, the 
Polywell tm might still function as a fusion reactor by using a larger radius R. 

(2) Initially monoenergetic electron or ion distributions may be driven toward 
thermodynamic equilibrium (approaching Maxwellian distributions) by enhanced 
collisional mechanisms, such as the beam-plasma and two-stream instabilities [110]. 
It would be desirable to prevent the accumulation of cold electrons, in order to 
avoid the beam-plasma instability. 

(3) The polyhedral set of point cusps has magnetic field nulls on its surface, 
which might not confine high-pressure plasma well. However, no significant particle 
losses were associated with the null points in the early operation of the spherically 
convergent ion focus (SCIF) experiment [111]. 

2.11. Modified Penning trap 

Barnes and Turner [112] found that a time- varying external electrostatic quadrupole 
field in a uniform static magnetic field could induce large-amplitude, nonlinear 
oscillations of a non-neutral plasma with very high plasma densities at the centre in 
one-, two- and three-dimensional cases. These oscillations are predicted to be stable 
at large amplitude. 

By definition, the Brillouin ratio is the ratio of relativistic electron energy density 
to magnetic field energy density: 

Brillouin ratio = 2 fi 0 nmc 2 /B 2 = n/n B (51) 

where n B is the Brillouin density. In an ordinary equilibrium case n/n B ^l, the 
‘Brillouin limit’. Turner and Barnes [113] found that a class of flows with deviatoric 
strain can support charged particle densities with Brillouin ratios locally exceeding 
unity. This configuration does not require large-amplitude plasma oscillations. They 
showed that, in a toroidal geometry with a purely poloidal magnetic field, local 
values of the Brillouin ratio n/n B > 2 could be sustained. The poloidal magnetic field 
could be produced by a toroidal quadrupolar current configuration outside the 
plasma, and a toroidal octupolar electric field applied outside the plasma could be 
used to sustain the required plasma configuration. With a 1/r term in the toroidal 
flow velocity, the Brillouin ratio at the x point (the centre of the plasma) becomes 
infinite. In contrast to the usual Penning-trap-confined non-neutral plasmas, this 
configuration is not in thermodynamic equilibrium. Large values of the Brillouin 
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ratio are feasible in a cold plasma when the electron velocity shear length 

/j s ■ n/ojp e , 

Barnes et al [114] studied non-neutral plasma confinement in a uniform magnetic 
field with an applied electrostatic quadrupole field (a Penning trap). They found a 
class of confinement that does not require a large velocity shear. The plasma 
reactivity is proportional to where ( ) denotes a space and time average. Since 
the Brillouin limit applies only in the volume- averaged sense, the reactivity can be 
enhanced by inducing large spatial variations of n over the plasma volume (such as 
ultrahigh n at the centre). Low-energy, low-P^ (canonical angular momentum) 
charged particles (such as electrons or deuterons) could be injected into this 
spherical Penning trap. The effective parabolic spherical potential well of the applied 
electromagnetic field would accelerate the particles radially inwards to energies 
—100 keV. Convergence near the centre (r = 0) produces a dense, inertially 
confined, non-neutral plasma core. Strong focusing occurs only for a single 
charge-to-mass ratio q/m , so a pure deuterium plasma might be used. Miyamoto et 
al [115] showed that, with a judicious choice of a rational integer determining the 
electric field shape, the DD fusion reaction products 3 H and 3 He could also be 
focused to pass through the central region, enhancing the fusion power density. 

Another possibility is an initially pure electron plasma, with deuterium-tritium 
ions trapped in the virtual cathode. In the electron-plasma version low-energy 
electrons would be introduced through one of the negative end-caps of the 
electrostatic quadrupole, as illustrated in figure 17. In a reference frame rotating at 
angular frequency Q = - eBjlm , the Poisson equation becomes [114] 

V 2 <S> = (efe 0 )(n n ~nl (52) 

The applied electric field is produced by electrodes with hyperbolic surfaces tangent 
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Figure 17, Schematic diagram of the electrostatic quadrupole Penning Trap [114]. 
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to the spherical boundary of the potential well. Particles with energies <2V 0 /3 will 
be confined, where the applied voltage V 0 = n B ea 2 1 4s 0 . 

The cold injected electron beam will be reflected near the centre of the spherical 
well ( r = 0). If P e = 0, then the characteristic radius r c of the dense central core could 
be comparable to a local Debye length. With finite perpendicular injection 
temperature T±, the angular momentum conservation equation yields 

r c ~r 0 (T ± /W) m (53) 

where W is the electron injection energy and r 0 ~a/ 6 is the radius of a virtual anode 
potential peak in the rotating reference frame. The electron beam velocity v varies 
slowly, except near the inner and outer turning points. From the steady-state 
continuity equation with F^O, the density varies radially as n(r)~n(a)a 2 /r 2 , 
where a is the radius of spherical bounding surface, and the central core density 
n c ** n B a 2 jr\. Since n 2 decreases rapidly for r > r c , the average of n 2 is estimated to 
be 

{n 2 ) ~ n 2 Q rl/a 3 ~ 16£oV 2 0 / e 2 a 2 r c . (54) 

This proportionality to V 0 /a 3 r c shows that the reactivity can be enhanced by high 
voltages and small radii. 

The scaling was confirmed by numerically computing self-consistent Vlasov 
equilibria for an example case with a = 3 mm, V 0 = 120 kv, electron beam radial 
injection at W = 50 eV, 7 = 0.125 A, T ± = 10 -4 eV and T r = 5eV. The numerical 
result showed that a density n c = 4 X 10 2S m -3 could be confined within r Q ~ 0.4 /xm 
[114]. For comparison, the above approximate equations yield r c = 0.71 pm, 
n R ** 3 X 10 18 m -3 , n c ** n B a 2 /r 2 = 5.3 X 10 25 m -3 and (n 2 ) « 3.7 X 10 40 m -3 . Thus, for 
this example case the approximate scalings yield a core radius and density not far 
from the numerically computed values. 

A second example case considered the effects of ions that partially neutralize the 
electron space charge. Assuming n-Jn -0.5 for r< 1 ju,m, n-Jn = (1-10 -6 ) for 
1 pm<r<r Q , n-Jn = 0 for r>r 0 and 7 = 0.425 A, similar potential and density 
profiles were obtained, with r c ~ 0.3 mm and « c «4x lO 25 m -3 . The neutral back- 
ground penetrates most of the plasma, but within a small ‘burn-out’ radius high ion 
temperatures may be sustained. The central ion confinement time is about 0.25 ms, 
and a neutral gas pressure about 1-10 mPa would be appropriate. For this example 
case the predicted central ion temperature was 7.8 keV, and the DT fusion power 
would be about 50 mW, and the energy gain Q ~ 10~ 3 . Calculations indicate that this 
configuration is stable against the two-stream instability, that Coulomb collisions do 
not significantly alter the electron beam velocity distribution, and that small angular 
momenta are attainable [114]. 

Although such a Penning trap would not produce a large total power output, it 
could produce useful fusion reaction rates in a laboratory scale device. Higher-order 
multipole traps could generate higher power outputs [114, 116]. 

Tiouririne et al [117] studied the use of higher-order multipole traps for 
non-neutral plasma confinement. The total electrostatic potential 4> is the sum of the 
externally applied potential <f> ap p and the potential <£ s generated by the plasma space 
charge. The plasma density, computed in terms of using the Poisson equation, was 
assumed to be zero outside a closed boundary, and the plasma source potential <E> S 
was computed assuming that it vanishes at large distances. Then the required value 
of applied potential # app = O - <J> S . The magnetic field contours, plasma density 



1568 


TJ Dolan 


Contours 

I = 2 (magnetic quadruple) 1 = 3 (magnetic hexapole) 

(a) isodensity (d) isodensity 
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electrostatic potential electrostatic potential 


z 


Figure 18. Contours of plasma density, magnetic flux, and externally applied electro- 
static potential for a magnetic quadrupole (a,b,c) and for magnetic hexapole { d,e,f ) 
1 117 ]- 

contours and externally applied potential contours computed for the magnetic 
quadrupole (l = 2) and magnetic hexapole (l = 3) cases in cylindrical geometry are 
shown in figure 18. The global Brillouin ratio was found to be £ = 21/(1 + 1), which 
approaches 2 for high-order multipoles. 
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Turner et al [118] calculated the frequency of the non-neutral plasma anharmonic 
oscillations and found that <al<a c varies from 2" 1/2 when S = 0 to (2/3) 1/2 when 8 = 1, 
where S = (p 2 ~ Pi)/(,P2 + Pi)> and Pi> P2 are the turning points of the oscillations of 
the plasma boundary. 

Although the studies of modified Penning traps are just beginning, they have 
shown several interesting features: 

(i) In a non-neutral plasma, densities exceeding the Brillouin limit can be 
confined by several means: in an oscillating plasma; in a flow with deviatoric strain 
(related to velocity shear) such as a torodial magnetic quadrupole with applied 
electric field; in a homogeneous magnetic field with an applied electrostatic 
quadrupole and low-energy particle injection; in a higher-order magnetic multipole 
field with an applied electrostatic multipole. 

(ii) A large degree of focusing is possible, with resultant very high densities n e 
inside the core radius r c . This focusing can result in high values of reactivity, which is 
proportional to ( n (i) 2 ). 

(iii) Higher-order magnetic multipoles can have global Brillouin ratios ap- 
proaching 1, 

While these results are not immediately conducive to an attractive fusion power 
reactor, further developments along these lines could result in newer configurations 
with more attractive properties. The evolution towards higher-order magnetic 
multipoles is similar to that of the electromagnetic trap experiments, which started 
with spindle cusps, then went to the toroidal quadrupole and higher-order 
multipoles, such as a linear set of ring cusps. 


3. MEPC experiments 

3.1. Spindle cusps and mirrors 

Early experimental studies of electrostatic plugging were conducted by Lavrent’ev 
[5, 36] and by Hilton et al [119]. Due to financial limitations, most early experiments 
were pulsed spindle cusps, although spindle cusps are poor configurations for 
electrostatic plugging. Electrons were injected from electron guns in the point cusps, 
with plasma production by ionization of neutral gas. Pulsed magnetic fields have 
several disadvantages for MEPC experiments: 

(i) eddy currents cause magnetic field errors; 

(ii) during magnetic field decay the plasma drift is outwards; 

(iii) plasma density build-up and heating may take longer than the magnetic field 
lasts [120]. 

Experiments with spindle cusp electromagnetic traps like S-l, S-3M, S-4M, and 
Jupiter-1 A showed that: 

(i) The electron confinement time increased by three orders of magnitude when 
the plugging voltage was applied. 

(ii) The electron confinement time was inversely proportional to the neutral gas 
pressure [12]. 

(iii) The rate of electron loss across the magnetic field was similar to that 
expected from classical diffusion due to collisions with ions and neutral atoms. 
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W, keV 


Figure 19. Measured electron energy distribution in the anode region immediately after 
cut-off of electron beam injection, and 1.5 ms later [81]. 


(iv) The plasma density increased linearly with applied voltage, up to a stability 
limit. 

(v) the plasma density increased proportional to B 7 [122, 123]. 

(vi) The mean ion energy scaled linearly with applied voltage [124]. 

(vii) The electron distribution function [81] shown in figure 19 is comparable to 
the theoretical curve of figure 15. 

(viii) A deep, negative potential well could be sustained for many ms, provided 
that the neutral gas pressure was low. 

(ix) Electron orbits and loss , angles agree with theoretically predicted values, 
provided that the starting point is within an adiabatic region [125]. 

(x) Diocotron oscillations occurred, as predicted theoretically. Experimental 
observations of diocotron oscillations in the ring cusp indicated the following: 

(1) The frequency of the oscillations was proportional to applied voltage and 
inversely proportional to magnetic field strength, as expected. 

(2) The growth rate of the instability was proportional to neutral gas pressure, 
and thus to the density of cold electrons produced by ionization. 

(3) By segmenting the ring cusp cathode, the relative ion currents as a 
function of azimuthal angle were measured. These measurements indicated that 
the long-wavelength instability had a mode number m — 1 or two wave periods 
around the circumference of the ring. Higher-mode numbers were apparently 
more stable [126]. 

(4) The anode region electron density remained below the short-wavelength 
diocotron limit of equation (27). 

Sugisaki [127] showed that the application of electrostatic fields to a magnetically 
confined plasma prolonged the plasma decay time. Strijland [128] studied particle 
orbits in an electrostatically plugged spindle cusp and the variation of confinement 
time with diaphragm diameter. With electrostatic plugging of a magnetic mirror, 
Nishida et al [129] observed a reduction of ion losses to 15% of the unplugged value, 
and ion heating by E X B rotation. 

Vdovin and co-workers [130] studied electron loss processes from the S-4M 
electromagnetic trap. They derived an equation for equilibrium plasma density that 
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varied with electron beam current, electron energy and neutral gas pressure in 
accord with experimental data. Plasma electron diffusion was dominated by 
scattering off gas molecules and molecular ions. They suggested that anomalous 
electron transport in the ring cusp region may be caused by a long-wavelength 
two-stream instability. The injected electrons had an anomalously fast loss rate in 
the S-4M electromagnetic trap, and neutral gas effects kept the plasma cold (a few 
eV temperature). 

According to measurements by Komarov et al [131] with a two-wavelength 
microwave interferometer in the Jupiter-IA spindle cusp device, the central density 
was an order of magnitude higher than the average density. Zalesskii et al [132] 
found that the the plasma accumulation rate with microwave heating was much 
faster than with electron beam injection, and about 50% of the microwave power 
was coupled to the plasma. 

Stepanenko and Komarov [133] measured the Jupiter-1 A electron and ion 
energy spectra during microwave heating with 10% power absorption, and deter- 
mined the electron energy confinement time to be about 0.9 ms. The dominant ion 
energy loss mechanism was charge exchange. 

Electron and ion energy spectra were measured with gridded electrostatic energy 
analysers at the point cusps. The electron energy distribution is shown in figure 19. 
The ion energy spectra were measured with a drift mass spectrometer combined 
with an electrostatic ion-energy analyser. The distribution appeared to have two 
Maxwellian components. The higher-temperature component was believed to be 
caused by the initial formation of the potential well. The average energy of the ions 
increased linearly with applied voltage, {Wj) ~ 0.18<j& A [36]. 

Komarov and Stepanenko [134] measured the spatial variation of the total 
number of electrons N c leaving the ring cusp gap, and found a narrow distribution 
with width S~p e , as shown in figure 20. At an electron energy of 0.55 keV 
p e = 0.1 mm. The distribution was narrower at high plasma densities, indicating that 
electron losses along the magnetic field became important relative to cross-field 
diffusion losses. The model profile with w = p e is roughly consistent with the 
experimental data, but w = 2 p e is used here in estimating future device parameters 
to avoid underestimating A<}>. Komarov and co-workers [135] also measured the 
distributions of electrons and ions in a point cusp gap, finding a narrow electron 
stream with radius <1 mm. 



2, mm 


Figure 20. Experimentally measured electron density distribution (full curve) in the ring 
cusp of the Jupiter-1 A device at $ A = 6 kV 5 = 0.82 T, and model profiles n/n a = 
1/[1 + (z/w) 2 ] with iv = p c (broken curve) and w = 2 p e (dotted curve). 
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In the KEMP-II spindle cusp experiment Stansfield and colleagues [136] 
measured n— 4xl0 17 m -3 , and electron energies — IkeV. Ionization of gas pro- 
longed the plasma density persistence after electron beam cut-off. 


3.2. Jupiter-1 M 

The Jupiter-IM spindle cusp experiment at Kharkov [137] had pulsed magnetic 
point cusp fields =£13 T (25 ms half-period), ring cusp fields =£4 T, electron injection 
voltages ~4 kV, plugging voltages <f> A — 6 kV, electron beam injection current 
£2mA, neutral gass pressure 0.1-1 mPa and plasma volume 30 cm 3 . It attained 
plasma densities up to 10 18 m -3 , with temperatures —100 eV and confinement times 
—1 ms. The short duration of the pulsed magnetic fields did not allow enough time 
for build-up of higher plasma densities. The plasma losses from the point cusps were 
less than from the ring cusp [138]. 

Komarov and colleagues [139] used gridded analysers to measure the energy 
distributions of electrons and ions emerging from the cusp gaps after removal of the 
plugging voltage. The mean ion energy observed under those conditions was about 
one-fifth of the plugging voltage. For a few microseconds after electrostatic plugging 
was removed, the remaining ions were accelerated to higher energies by the plasma 
potential [140]. The measured potential barrier in the ring cusp for ions — 32] 
exceeded the calculated barrier —27]. The barrier height in the point cusp was 
practically zero, but most of the ion losses occurred through the ring cusp. Point 
cusp ion losses were probably retarted by the centrifugal and magnetic field gradient 
forces [141]. 

Azovskii and colleagues [142] found that the inner bounding surface of the 
plasma in the Jupiter- 1M experiment has a radius r 0 = (2 Rp 0 ) m in the point cusp 
(where p 0 is the electron Larmor radius in the point cusp), as predicted from flux 
mapping to the ring cusp. The density of plasma flowing out of the point cusp was 
two orders of magnitude lower than the central plasma density, and A <f> ~ 0.5 kV 
was inferred from the shape of the ion energy distribution emerging from the point 
cusp [143]. 

Azovskii and co-workers [144] measured the plasma density in Jupiter-IM 
simultaneously by a microwave interferometer and by the integral of the ion loss 
current, and found that n ~ 10 ls mT 3 . The lifetimes of electrons injected from the 
gun were much shorter than the lifetimes of plasma electrons produced by 
ionization. The confinement time of the bulk plasma was about a factor of five lower 
than classical, assuming that the electrons were produced by ionization. The injected 
fast electrons left quickly (—0.01 ms), perhaps due to a two-stream instability, and 
imparted only about 1-2% of their energy to the plasma. The electron current had 
0.1-1 MHz modulations. Probe measurements showed that the electron density 
profile was not consistent with a theoretical probfile based on classical transport 
[145]. 

Factors influencing the electron confinement times included the magnetic field 
geometry, impurities in the plasma, ionization of neutral gas and electron mobility in 
the electric field. The effect of the centrifugal force on radial electron transport 
would be less significant in large devices. Since the electron mobility is related to the 
neutral gas density, its contribution was significant in small devices, but would be 
reduced in high-temperature plasmas, where most of the neutral gas becomes 
ionized [146]. 
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From ion energy distributions measured outside the point cusp and ring cusps of 
the Jupiter-IM device at <£ A = 6kV, Azovskii and co-workers [147] deduced the 
amount of potential sagging in the anode gaps, and found A<f> (ring cusp) = 0.4 kV, 
A cf> (point cusp) = 0.6 kV. These values correspond to uncompensated anode 
electron densities n a ~ 10 17 m -3 , while the central electron density n ~ 10 18 m -3 . The 
electronic space charge in the anode gaps was partially compensated by ions flowing 
out the cusps. Although the ion potential barrier fa was lower in the point cusp, the 
ion loss rate there was lower than in the ring cusp. From gridded analyser data they 
found that the ion energy distribution was nearly Maxwellian with 7j = 40-60 eV. 7] 
was slightly higher at the plasma edge, where the electric field was large [148]. 

Azovskii and colleagues [149] measured the electrostatic potential and electron 
temperature in Jupiter-IM. At the end of electron beam injection the potential sag 
in the point cusp cathode A<£ ~ 1 kV. The electron temperature reached 25 eV, then 
dropped rapidly after electron beam injection cut-off. Using an 8 keV Ar + probe 
beam, they measured the energy spectra of Ar 2+ ions emerging from the ring cusp 
and deduced the time variation of the central plasma potential. When energetic 
electrons were injected into low-pressure hydrogen, the plasma potential (fa + A 4>) 
initially was very large. As ionization built up the plasma density, the plasma 
potential decreased to a plateau value on the order of 0.1 $ A . After the electron 
beam injection was turned off, the plasma density and potential decayed to zero in a 
few milliseconds [150]. 

3.3. Magnetic multiples 

Nakamura et al [151] observed enhancement of ion confinement in a multipole 
magnetic field by a negative electrostatic potential well produced by electron 
injection into the plasma. Hershkowitz et al [152] observed similar effects. 

Lavrent’ev and colleagues [153, 154] studied the plasma produced by injection of 
100 mA of 1 keV electrons into a pulsed multipole trap. They measured n ~ 
10 17 m -3 , confinement time = 50 p,s and 7] =* 300 eV. Using the attenuation of a 
calcium ion beam during a 30 ps period of electron injection, they observed, a rapid 
density build-up to 4 x 10 17 m -3 , followed by a decay after electron beam cut-off. 
With probe measurements they found a uniform central plasma density dropping off 
steeply in the boundary magnetic field layer [155]. Lavrent’ev et al [156] found that 
the plasma lifetime in a multipole electromagnetic trap depended very strongly on 
the magnetic field in the gaps, decreasing from 1.3 ms to 0.1ms as the field was 
reduced from 0.3 T to 0.15 T. 

Moir et al [157] proposed a toroidal quadrupole MEPC experiment using soft 
iron to enhance the magnetic field shape and to minimize the required coil currents. 

Maffei et al [158] constructed a spherically convergent ion focus (SCIF) device 
based on Bussard’s PoIywell tm concept (a polyhedral set of point cusps), with 
planned values of B = 0.2 T, <p = 20 kV, I = 75 A from three electron guns and 
plasma production by electron cyclotron heating. 

3.4. ATOLL 

The ATOLL toroidal quadrupole magnetic cusp device is shown in figure 21. 
Typical parameter values are 2? a = 0.8 T, <i> A ~ 2 kV, p = 0.2-0.5 mPa hydrogen, 
impurity fraction = 1-2% (dominant mass = 28), n = 3 X 10 17 m -3 , T e = 100 eV and 
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Figure 2L The ATOLL toroidal quadrupole cusp device [160]. 


average ion energy = 100 eV. Ioffe et al [159] checked the alignment of the magnetic 
field with the anodes by observing the flow of electrons from an electron gun placed 
in the plasma region. The alignment accuracy was better than 0.2 mm. Diagnostics 
included an 8 mm microwave interferometer, currents drawn by electrodes, a 
charged particle energy analyser, Langmuir probes (including a tiny probe sus- 
pended by a 12 fxi n tungsten fibre), a charge exchange neutral atom energy analyser, 
a pulsed vacuum guage and a partial pressure analyser. 

In one mode of operation the plasma was sustained by electron injection from a 
gun plus secondary emission at the reflecting cathodes. This mode resembled a 
Penning discharge with the end walls outside the cusps. 

In another mode of operation Ioffe and co-workers [160] minimized the plasma 
impurity content by careful wall preparation and by the use of very pure hydrogen. 
Ionization of neutral gas made up for particle losses, and the electrons were heated 
by the intense electric field in the boundary layer. At high magnetic fields and low 
pressure the discharge could become self-sustaining without electron injection or 
secondary emission from the cathodes. This regime had similar plasma parameters, 
but less noise, better ion confinement and higher % than the first regime. This 
‘emission-free’ discharge was estinguished by adding 5% ,pf an impurity, such as 
nitrogen. Impurities increased radiation losses and diffusion rates, and reduced the 
potential well depth. 

The plasma potential increased linearly with magnetic field and had a nearly 
parabolic radial distribution. Ions produced by ionization at the edge were heated as 
they fell into the potential well, and the average ion energy increased with increasing 
magnetic field. From the plasma confinement time it was inferred that the potential 
barrier impeding ion loss through the cusps was very small. Ioffe and colleagues 
[161] found that the electrostatic potential in the cusp gaps was equal to the 
potential at the plasma centre, to within experimental accuracy (10%). The 
hypothesis of a narrow ‘neutral corridor’ through the cusps with no confining 
potential barrier (<f > j = 0) was also consistent with the measured variation of plasma 
potential with magnetic field. 

The electron potential barrier <j> e ~ 107^, so electron losses were predominantly 
by cross-field diffusion. From the ionization rate Ioffe and co-workers [162, 163] 
inferred a plasma co nfi nement time of 0.18 ms. This was close to the theoretical 
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confinement time, based on the ion acoustic instability, of 0.12 ms, and much shorter 
than the classical confinement time of 50-150 ms. Similarly, the measured electron 
density profile was close to that predicted above for the ion acoustic instability, with 
a gradient scale length on the order of the ion sound Larmor radius, and far from 
the profile based on classical diffusion. Piterskii and co-workers [164, 165] observed 
oscillations at 40-80 MHz propagating in the electron diamagnetic drift direction at 
about 10 s ms -1 , which appeared to be consistent with ion acoustic instability theory. 
According to Lavrent’ev [166] the cross-field transport resulting from the observed 
profiles and oscillations in ATOLL should be inwards. 

Later measurements revealed that three instabilities were occurring simul- 
taneously in ATOLL: diocotron oscillations in the cusp gaps, a long-wavelength ion 
acoustic instability in the plasma interior and a lower hybrid drift instability at the 
plasma edge. 

Piterskii et al [167] observed diocotron oscillations in the cusp gaps at frequencies 
of 40-70 MHz, with wavelengths of 1-2 era and peak amplitudes -T00V. The wave 
phase velocity was somewhat lower than the electron E.XB drift velocity in the 
anode gaps, by a factor of 1.5-2; but this difference could be partly due to 
simplifying assumptions used in the analysis, such as the shapes of the electron 
density distribution and electric field distribution in the anode gaps. The theoretical 
maximum growth rate of the diocotron oscillations occurs at A ~ 16a, which is within 
the range of wavelengths observed. The instability became stronger with increasing 
plasma potential and decreasing magnetic field, as expected for diocotron oscilla- 
tions. In the central plasma the amplitude of the diocotron oscillations decayed with 
distance from the anode gaps with a scale length of about 5 cm, as expected 
theoretically. 

Piterskii and co-workers [168, 169] found that the lower hybrid drift instability is 
dominant in the boundary region where the density is low. The theoretical values of 
most unstable wavelength (1cm), frequency (12 MHz) and phase velocity (1.2 X 
10 s m s -1 ) are consistent with experimental data. This instability causes anomalous 
electron transport, flattening the electron density profile until the density gradient 
scale length (rc/Vn) exceeds the local ion Larmor radius pi- Thus, the confinement 
times in ATOLL were much shorter than classical confinement times, and the 
anomalous transport was well understood in terms of the lower hybrid drift and 
long-wavelength ion acoustic instabilities. 

3.5. Jupiter-2M 

The proposed Jupiter-2 experiment [170] was to be a linear set of ring cusps with 
B*=5T (point cusps), </> A = 25kV, V = 0.5m 3 , n>10 19 m -3 , 7j~lkeV and r E ~ 
100 ms. Due to financial limitations, a one-third scale model of Jupiter 2, called 
Jupiter-2M [171] has been built first (figure 22), with a primary goal of studying 
electron transport. 

The Jupiter-2M device has seven ring cusps with R = 0.215 m in the central 
region flanked by one small ring cusp and point cusp at each end. The axial length 
between the point cusps is 1.30 m. The coils are pulsed from a capacitor bank and 
reach peak fields in about 33 ms. The maximum attainable ring cusp and point cusp 
magnetic fields are about 1.94 arid 3.88 T, respectively, but lower fields have been 
used in experiments. The shape of the vacuum magnetic field is illustrated in figure 
23 [171]. 
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Figure 22. The Jupiter-2M device (left) and electrode region detail (right) [171]. 


vacuum pumping port 


The plugging electrode positions were adjusted by screws until comb probes at 
three azimuthal locations 120° apart showed no asymmetry of the electron loss 
currents emerging from the cusps. The total plasma volume is 50 1, and the volume 
of the uniform central region is 16 1. 

The classical electron diffusional loss rate is [172, 173] 

hm = [Ai(l + 0.5 <£ p / Teo) + D ei ]Fn 0 R 2 electrons/s (55) 

where £> a is the electron diffusion coefficient due to collisions with neutral atoms, 
d> p = </>; 4- Ac/), D ci is the electron diffusion coefficient due to collisions with ions, n 0 is 
the central electron density, T e0 is the central electron temperature, and F is the 
geometrical factor having dimensions of reciprocal length. The 4> p term represents 
electron mobility in the strong edge electric field. The electron confinement time is 

Tjiir — NJ i d i ff — l.HjzR 3 n 0 / /<j ^ (56) 

where N c is the total number of electrons in the plasma. For the geometry of the 
Jupiter-2M experiment F = 900/cm. 



Figure 23. Vacuum magnetic field lines in the Jupiter-2M device [171]. 
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n.. 10 1S m’ 3 


Figure 24. Ratio of experimental to theoretical electron loss currents, as a function of 
plasma density. B = 0.375 T (+), 0.5 T (•), 0.625 T (*) [174). 


Preliminary low-field operation at B ~ 0.5 T (ring cusp) with a 1.3 kV, 0.12 A 
electron beam injection into 2 X 10 -6 Torr (H 2 ) produced a plasma with n = 
5 X 10 17 m" 3 , T e = 13 eV, T e > 2 ms. With <£ A ~ 2 kV, the plasma potential was about 
—300 V relative to the grounded chamber walls. About 80% of the injected electron 
beam energy appeared in the plasma. Although the point cusp anode diameters 
were 25 mm, the ion loss through the point cusps was only 6% of the total ion loss 
rate. 

High-frequency oscillations were not observed in Jupiter-2M. The large uniform 
central plasma may serve as a buffer to damp out surface oscillations. Vdovin et al 
[174] found that the ratio of experimental electron loss current to theoretical loss 
current decreases with increasing plasma density, as shown in figure 24. At higher 
plasma densities, the experimental loss rate was about twice the classical value. The 
measured ion loss rate corresponds to an ion potential barrier 

Auxiliary electrodes were used to suppress secondary electron emission from the 
plugging cathodes. Diaphragms limiting the plasma boundary in the point cusps and 
between the ring cusps served two functions: (i) to reduce particle flow onto the 
anodes, where gas emission stimulated by particle bombardment would tend to fill 
the gap with cold plasma, and (ii) to study the influence of the diaphragm size on 
plasma transport. Ring diaphragms with diameters of 21, 24 and 28 cm corresponded 
to ring cusp diffusion layer widths a = 0.9, 1.6 and 2.0 mm. 

The theoretical electron density profile depends little on the temperature profile. 
The electron density profile is flat in the centre and falls nearly linearly to zero near 
the plasma edge. There is a uniform central plasma with radius —5 cm and length 
~80 cm. As the diffusion layer width a is increased from 0.9 to 2.0 mm, the electron 
confinement improves, but ion confinement degrades, due to the increase of A<j>. 
There is also an increased probability of arcing with large gap widths, due to anode 
bombardment and outgassing. At the optimum width a opi = 1.6 mm the maximum 
plasma density n = 10 18 m“ 3 is attained. The electron confinement time deduced 
from the total number of electrons in the plasma divided by the total electron 
diffusion loss current is z e = 3 ms. Future plans are to increase the magnetic field to 
1-2 T, and to inject hot plasma with a plasma gun. 
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3.6. Impurities 

Impurity effects in MEPC devices have been studied very little. Lavrent’ev [94] 
showed the/ theoretical effect of increasing Z efl on the energy confinement time. 
However, nobody has computed impurity source terms, transport and build-up. 

A continuous impurity accumulation was observed in the ATOLL device during 
the injection pulse. A typical impurity level was 1-2% with the dominant species at 
mass 28 [159]. The addition of a few per cent impurity extinguished the ‘emission 
free’ regime of ATOLL [160]. There have been almost no studies of impurity 
profiles, impurity transport, radiative power loss or methods of impurity removal 
from an MEPC plasma. Impurity studies need greater emphasis in the future. As 
Gormezano [39] stated: ‘Any electrostatic system will constitute a very good 
containment system for heavy impurities. The presence of electrodes in environ- 
ments with high electric fields and high particle fluxes may constitute a very large 
source of impurities. This might ultimately be the greatest problem’. 

3.7. Comparison of experiments 

Key parameters of various MEPC experiments are summarized in table 2. In some 
cases the effective temperature is taken to be two-thirds of the average particle 
energy. The early experiments at Kharkov (S-l, S-3, S-3M and S-4M), Quebec, 
Missouri and Jutphaas were small spindle cusp plasmas with weak fields, and neutral 
gas effects were dominant. Such experiments are of less interest now, but they 
served to demonstrate the principles of MEPC. The multiple point cusp SCIF device 
was shut down by budget cuts before reaching its projected operating parameters, 
but it may be restarted at another laboratory. 

The - Jupiter-1 A device achieved good ion and electron temperatures (0.4- 
1.5 keV). It also demonstrated the effectiveness of electron cyclotron heating, 
although at low powers. The Jupiter-IM experiment was not able to fully utilize its 
high magnetic fields, because the duration of the peak field was too short to permit 
accumulation of high plasma density. Due to their low ratios of plasma volume to 
surface area V/SR =£0.03, spindle cusps are not worthwhile considering for future 
experiments. 

The linear multipole devices demonstrated confinement in that geometry and 
validated theoretical predictions of critical loss angles. Due to end losses, these 
devices did not achieve high plasma pressures. 

The ATOLL totoridal quadrupole was carefully built and operated, but the 
electron energy transport was anomalous. Detailed studies related the anomalous 
transport to the lower hydrid drift instability (plasma edge) and to the ion acoustic 
instability (plasma interior). ATOLL was shut down in the 1980s. 

The Jupiter-2M experiment has attained the most promising results. Its demon- 
stration of near-classical transport is stimulating interest in MEPC research. 

3.7.1. Density , temperature and confinement time scalings. The plasma density in 
spindle cusps was typically 10 17 -10 18 m -3 ' The attainment of higher plasma density 
has been a major concern for the reactor prospects of MEPC devices. The plasma 
density appears to scale proportional to B 2 [122], provided that the magnetic field 
duration is long enough for density accumulation. The central plasma density n is 
found to be significantly higher than the electron density n a in the anode region, but 



Table 2. Parameters of some MEPC experiments. Italicized values are estimates* 
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more definitive measurements of the ratio njn are needed at high plasma 
temperatures. 

The measured electron energy distributions resemble shifted Maxwellians, as 
expected. Computer models indicate that, if neutral gas effects are not dominant, the 
plasma temperatures scale proportional to the applied voltage. Measured ion 
temperatures have scaled linearly with 4> A . For the plasma density and temperature 
in MEPC devices, we assume the following scaling 

n<xB 2 (57) 

r e ~7]~0.054> A . (58) 

The scalings were used to estimate the hypothetical reactor parameters in tables 1 
and 2. That the temperatures in Jupiter-2M were below this assumed scaling may be 
due to the short duration of the peak magnetic field, the weak heating power and 
neutral gas effects. 

For small spindle cusp devices the confinement times have been T~ 0.1-1 ms. 
Little distinction has been made experimentally between particle confinement and 
energy confinement. For ring cusps the confinement time is assumed to scale roughly 
in accordance with equation (40) 

r d j ff =0.09T ei ai?/p b p a . (40) 

This equation predicts r diff ~ O.S ms for Jupiter-2M parameters of table 2. The 

experimentally measured confinement time was 2-3 ms. (If the central electron 
temperature were higher than the 15 eV reported from a Langmuir probe under the 
magnet coils, then the theroretical T diff would be longer.) A more precise 

computation specific to the Jupiter-2M geometry [173] predicts classical confinement 
times about twice the measured value. 

3.7.2. Dimensionless parameters. Three dimensionless parameter ratios are listed in 
table 2: alp,., /3 C , and v a -r/R. The a!p e represents the number of electron Larmor 
radii within one-half of the ring cusp anode gap. In MEPC experiments this number 
has ranged from 5 to 100, and for a reactor it would probably be in the order of 
20-30. The optimum anode gap width is of the order of 1.5-2 mm, and the electron 
Larmor radius is fixed by the magnetic field and electron temperature in a reactor to 
values of the order of 30-100 p. m. 

The second dimensionless parameter is the ratio of average plasma pressure to 
the peak magnetic field pressure at the coil 

B e ^2p 0 nk(T e + TI)lBl (59) 

where B c is the peak magnetic field at the coil ~ B a . (In magnetic confinement fusion 
it is more common to define beta relative to the vacuum magnetic field at the plasma 
centre, but that central magnetic field is zero in magnetic cusps, which would yield 
an infinite value of beta.) The values of /3 C in past MEPC experiments have been 
very low, of the order of 0.001 to 0.1%. Using the density and temperature scalings 
of equations (57) and (58) results in [3 C = 4.5% for the hypothetical reactor of table 
2. For comparison, the estimated value for the ITER tokamak (Troyon coefficient 
g = 3, / = 24 MA, minor radius = 3m, B =6T, B C = 13T) is /3 C = 0,85%. 

The third dimensionless parameter in table 2 is v e r/R , which is the ratio of the 
confinement time to the approximate electron transit time through the plasma. 
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Present experiments attain confinement values of about 10 4 transit times. Without 
effective plugging, the plasma would only be confined for tens of transit times. 

This number of transit times is relevant to a mode of free plasma flow out cusp 
gaps. It is less relevant to MEPC devices, because ions are confined electrostatically, 
and their confinement times are roughly equal to the electron cross-field diffusion 
time. Nevertheless, it is instructive to show that, if the ions were freely flowing out 
cusp gaps, those gaps would be very narrow. According to the free-flow model, the 
ion confinement time is 

z~r 2 p LjuR8. (6) 

For Jupiter-2M, R = 0.22 m, r p ~ 0.1 m, L x = 0.105 m, u = 7.6 X 10 4 ms _1 , and t * 2 
3 ms. Using this equation we estimate the effective cusp gap half- width 5 ~ 
2.1 x 10 -5 m. Since the electron Larmor radius p a = 2.2 x 1CT 5 m, it is seen that 
8 » p a . If free ion flow were occurring, it would be through a very narrow gap. 

3.7.3. Anomalous transport. The diocotron instability in the ring cusp regions is 
believed to play a benficial role by removing anode-trapped cold electrons, which 
tend to screen out the anode potential, but other instabilities have deleterious 
effects. The ATOLL toroidal quadrupole cusp experiemnt had anomalous transport 
(hundreds of times faster than classical) in all operating regimes, associated with the 
long-wavelength ion acoustic instability and the lower hybrid drift instability. The 
Jupiter-2M multiple ring cusp experiment, however, has demonstrated nearly 
classical electron transport. So far, the reason for the difference is not understood. 
There are several hypotheses, as follows: 

(i) In Jupiter-2M directed electron injection from thermoelectric cathodes was 
used. According to Lavrent’ev [1], electron injection in ATOLL by secondary 
emission from the cathodes introduced electrons into the diffusional layer (rather 
than into the central Btillouin region), which could foster the two-stream instability. 
The rate of electron production was two orders of magnitude greater than the 
electron loss rate by classical diffusion, so classical diffusion could not accommodate 
the required outflow rate of spent electrons from the trap. 

(ii) ATOLL has two-dimensional curvature and only toroidal symmetry in the 
anode regions, while Jupiter-2M has one-dimensional curvature and (since h»a) 
approximately two degrees of symmetry. This lack of symmetry in ATOLL might be 
more conducive to growth of the harmful instabilities, such as the lower hybrid drift 
waves; or it might somehow tend to. reduce the beneficial effects of diocotron 
oscillations on cold electrons within the anodes. 

(iii) The large volume of field-free plasma in Jupiter-2M may serve to damp the 
growth of instabilities. 

More study will be required to falsify or confirm these hypotheses. 

3.7.4. Assessment. MEPC experiments typically have had small radii R ~ 10-20 cm 
and plasma volume <11, except for Jupiter-2M (50 1) and the SCIF device. The 
anode gap half-width ranges from 1 to 2 mm. A larger gap could permit excessive 
voltage sag A <j>, and a smaller gap would decrease the electron diffusion loss time. 
The magnetic fields have generally been below IT, except for a few pulsed 
experiments. The applied voltage has been below 10 kV. 

Many excellent measurements have been performed in MEPC devices. Detailed 
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probe measurements have shown that the electron density profile in the ring cusp 
anode has a half-width about one Larmor radius. The electrostatic potential profile 
has been mapped out with pulsed ion beam transmission measurements and gridded 
energy analysers. A deep negative potential well has been sustained for many 
milliseconds. The electron velocity distribution resembles a shifted Maxwellian, as 
expected. These measurements verify the basic theory of MEPC. 

The scalings of density, temperature and confinement time have been established 
within the current ranges of parameters. The confinement in some experiments is 
nearly classical, when neutral gas effects are taken into account. The diocotron 
oscillations in the ring cusps have been thoroughly studied and are well understood. 
The plasma waves associated with anomalous transport in ATOLL have been 
correlated with theoretical predictions. 

Why some devices have anomalous transport and others do not is not yet 
understood. More measurements are needed of the density and temperature 
profiles, impurity distribution and transport, neutral gas distribution, radiation losses 
and heat propagation. Two-dimensional computer simulations are needed to provide 
a link between analytic theory and experimental data. 

A common thread of deficiencies links almost all the past MEPC experiments: 
either weak or short-duration magnetic fields; relatively low voltages (<10kV); 
weak plasma heating; incomplete diagnostics; inadequate computer modelling. 
These deficiencies have been caused by low budgets, not by failures of the 
experimenters. 

Table 3 compares the parameters of Jupiter-2M with the proposed Jupiter-2 and 
Jupiter-T devices. If successful, Jupiter-T could demonstrate the attainment of 
thermonuclear parameters in a modest-sized device. 


4. Reactor concepts 

4.1. Reactor design studies 

Various cusp geometries suitable for electrostatic plugging were illustrated in figure 
8. Mirror geometries, such as a yin-yang field, are also possible for electrostatic 
plugging, but the plasma volume inside the anode-bounded flux surface is smaller 
than for cusp systems. 


Table 3. comparison of Jupiter-2M, Jupiter-2, and Jupiter-T experiments [174]. 


Parameter 

Jupiter-2M 

Jupiter-2 

Jupiter-T 

Ring cusp radius, R (m) 

0.215 

0.67 

1.4 

Length between ends, L (m) 

1.4 

3.9 

6.4 

Ring cusp magnetic field, £ a (T) 

0.5 

2.5 

4.0 

Point cusp magnetic field, B z (T) 

1 

5 

8 

Applied voltage, <j> A (kV) 

2.5 

25 

100 

Electron injection current, I c (A) 

0.2 

5 

30 

Plasma density, n (m"" 3 ) 

0.8 X 10 1S 

3 X 10 19 

Itf 20 

Ion temperature, T (keV) 

0.03 

1 

10 

Plasma volume, V (m 3 ) 

0.05 

0.5 

4 

Electron diffusion current, / d{ff (A) 

1.6 

14 

120 

Electron lifetime, c (s) 

0.002 

0.36 

1.3 
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Laurent’ev [175] considered several reactor geometries for MEPC reactors: a 
sphere, a cylinder and a torus. They would typically operate at 5 a = 7T, <j> 2 = 
200 kV, n ~ 10 20 m -3 , 71 — 30 keV, V — 800 m 3 and power density —4 MW m -3 . They 
would provide intense sources of neutrons, which could be used for production of 
tritum or plutonium. 

Pastukhov and Yushmanov [176, 177] estimated the power gain of an MEPC 
reactor. At each value of applied potential <p A there is a plasma temperature that 
ma ximiz es Q. The optimum gap width makes Sd>j4> A ~ 0.25. When the gap is wider, 
&<f> becomes large, reducing the plugging potentials <j> e and 4> r When the gap is 
narrower, electron diffusion losses are too great. At the optimum gap width, the 
ratio of fusion power to input power for a ring cusp plasma reduces to 

Q ~ 10~ 5 cf> 2 A B a V/S(l + nJn A ) m (60) 

where </> A is in kV and n t is the trapped electron density in the anode region. 
Assuming njn a ~ 0.3 and V /SR ~ 0.3 in the example case of table 1, it is found that 
Q — 10. This analysis did not include energy losses due to neutral gas and impurities, 
or the beneficial effect of the alpha-particle emission current helping to reduce the 
required cathode current. Alternatively, using the parameters of table I with Z eff = 3 
in equations (4D1) and (3F14) of [44], we find Q — 13. 

A spherical cusp, shown in figure 8, would be advantageous from a physics 
standpoint [178, 179], The sphere has the best surface-to-volume ratio, and the 
spherical electrode arrangement might provide some spherical focusing of charged 
particles. If effective, such focusing could increase the central plasma density and 
decrease the required magnetic field in the cusps. However, particle scattering by 
waves and by Coulomb collisions and reflections off the scalloped boundary would 
tend to diminish the degree of focusing. It would be difficult to align the electrodes 
with the magnetic field lines in a spherical reactor. 

A torodial hexapole reactor design had B a = 8 T, 4> A — 300 kV, n = 10 20 m -3 , 
T = 17 keV and a neutron wall loading = 1.6 MW m -2 [84]. The direct capital costs 
estimated for this reactor were comparable to those for tokamak and tandem mirror 
reactors [180]. 

From the results with Jupiter-2M, it appears that ion losses from the point cusps 
are not ai serious problem, so a linear set of ring cusps would be feasible for a 
reactor. Such a reactor would be superior to either the spherical cusp or toroidal 
hexapole cusp, because of its simplicity, lower structural mass and easier main- 
tenance. A linear MEPC reactor, as illustrated in figure 25, could have the 
parameters of table 1, with Q « 10, DT fusion power density = 3 MW m -3 , neutron 



Figure 25. One end of a linear MEPC reactor. 
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Figure 26 . Radiation collimator and voltage divider for high-voltage electrodes. 

wall loading = 3.6 MW m -2 , L = 50 m, and total fusion power = 4.2 GW. Radiation 
collimators and voltage dividers would be used to prevent high voltage breakdown, 
as illustrated in figure 26. 


4.2. Applications 

Besides power production, there are other potential applications for MEPC devices. 
Their good confinement of high-Z ions would make them well suited for stripping 
and confining heavy ions for accelerators, as suggested by Stix [181, 182]. 

Nadler et al [183] have suggested that small inertial-electrostatic confinement 
(IEC) devices operating in the glow discharge mode at 30-50 kV could become 
useful low-cost neutron sources producing >10 6 DD neutrons/s for a variety of 
applications, including detector calibration, diagnostics and training. By increasing 
the voltage, current and size, and using DT fuel, yields of 10 10 neutrons/s might be 
attainable. Such neutron sources could be stronger than 2S2 Cf sources (2.6 year 
half-life), and they would be simpler and cheaper than accelerator-based neutron 
sources. 

Miley et al [184] found that IEC device operation in the ‘star’ discharge mode 
maximizes the effective grid transmission factor for ions. An IEC device using 
D- 3 He fuel could generate 14.7 MeV protons for production of positron-emitting 
radioisotropes that are useful for positron emission tomography facilities. Miley et al 
[185] proposed the use of high-voltage collectors in a ‘Venetian blind’ configuration 
to convert the D- 3 He reaction product kinetic energy directly into electricity. 

Miley et al [186] studied hydrogen production by radiolysis of water with 
14.7 MeV protons from an IEC device. The 3.6 MeV alpha energy would be 
converted directly into electricity. The protons would pass through a thin chamber 
wall and interact with steam flowing through the cell, splitting about 13 water 
molecules per 100 eV of energy deposited by the protons. The steam would pass 
through a turbine to generate electricity, and then the hydrogen would be removed 
from the condensate. The overall system efficiency could be about 50%. 

Bussard [187] suggested the use of an electron-beam-heated propellant thruster 
powered by an electrostatic plasma device for space propulsion. Miley et al [188] 
suggested the use of IEC engines with Venetian blind energy convertors to power 
five 68-N arcjet/magnetoplasmadynamic thrusters with a specific impulse of 3000 s. 
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In pulsed operation, an IEC reactor with R = 5.25 m and mass = 5 1 could deliver 
12 MW of power to the thruster [189]. Bussard [190] proposed the use of a Polywell 
reactor for space propulsion. 

An MEPC plasma with T = 10-20 keV would be a good target for neutral beam 
injection at IV' = 150 keV with a beam ion fusion probability ~ 10 -2 [191]. Such a 
device could serve as a powerful neutron source for isotope production, medical 
therapy or materials testing. 


5. Issues of MEPC 

Although some recent confinement results appear promising, there are still many 
issues to be resolved before MEPC devices can become fusion reactors, as discussed 
below. 

5.1. Self-shielding of applied potentials 

The potential drop in point cusp anodes is very large, but point cusp ion losses can 
be kept low by the centrifugal and magnetic field gradient forces, as demonstrated in 
Jupiter-2M. At B *£ 6 T, values of n a =£ 1.5 X 10 19 m -3 are required to keep A<f> 
100 kV in the ring cusps. The potential sagging in the anode means that high applied 
voltages are needed to ensure substantial and <£;. 

5.2. Neutral gas effects 

Neural gas increases cross field transport rates and promotes high-voltage break- 
down. The plasma tends to act as a vacuum pump, and in future large MEPC 
devices the neutral gas effects should be less dominant than in small spindle cusps. 
Nevertheless, there will always be substantial recycling of hydrogen isotopes and 
helium off the walls, so shields should be used to prevent excessive flow of neutral 
gas to the anode regions, where build-up of cold trapped electrons would increase 
A<f>. High-speed vacuum pumping will be required in the electrode region. 


5.3. Electrodes 

Optimized plasma confinement by careful alginment of each anode was demon- 
strated in Jupiter-2M. Radiation shields and voltage dividers will be needed to 
sustain voltages <5 A >200kV. To avoid warping, the anodes should be actively 
cooled and made from materials with low thermal expansion. Although high-voltage 
holding in radiation fields must still be demonstrated in MEPC devices, the PSP-2 
rotating plasma experiment [192] has demonstrated the sustainment of 400 kV 
applied voltage with hot plasma present. 


5.4. Transport 

The reason why Jupiter-2M had good confinement and ATOLL had anomalous 
transport is not yet understood. Transport in MEPC devices has been studied much 
less than transport in tokamaks and tandem mirrors, so this lack of understanding is 
not surprising. Several hypotheses have been suggested: the effect of the electron 
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injection method, the effects of curvature and symmetry, and the effect of a large 
volume of field-free plasma. The extension of the good confinement to larger devices 
remains to be tested experimentally. A more powerful experiment is being designed 
to test the scaling relations at higher values of plasma density, temperature and 
confinement time. 


5.5. Power gain ratio 

The estimated power gain ratio for an MEPC reactor with nearly classical fuel ion 
confinement (but no alpha-particle confinement), <f> A = 400 kV, £ a = 6T, R = 4 m, is 
Q ~ 10. The power gain ratio could be increased by the following means: increasing 
the applied voltage; increasing the magnetic field; increasing the plasma volume; 
improving alpha energy co nfi nement. 

One means of improving alpha co nfin ement would be to add a modest toroidal 
field to a toroidal multipole device (similar to the Tormac configuration [193]). The 
toroidal field would help to confine alpha-particle energy, facilitating ignition, and 
the resultant magnetic shear might help to stabilize the lower hybrid drift instability 
in the plasma edge layer [194]. Another possibility might be the addition of a weak 
longitudinal magnetic field underneath the surface cusp field of a linear set of ring 
cusps. 


5.6. Plasma impurities 

Since alpha particles are poorly confined, helium ash accumulation could be slower 
than in tokamaks. However, sputtered wall atoms which penetrate into the plasma 
core would be well confined. Impurity ion accumulation would dilute the fuel 
density, cool the electrons by radiation losses and diminish the potential well depth. 
One means for impurity removal might be to use the centrifugal force of their 
rotation in a weak longitudinal field to expand their orbits so that a pumped limiter 
could extract them. Another possibility would be to expand the point cusp hole at 
one end of a linear reactor and to flow in fresh fuel from the other end, so that the 
weak pressure gradient gradually moved impurities to the hole. If effective impurity 
removal methods could not be found, then periodic shutdowns might be required for 
impurity removal, followed by restarts with fresh fuel. 


5.7. Materials issues 

The wall activation, afterheat, accidental release of radioactivity and waste disposal 
issues of MEPC reactors would be comparable to those of tokamaks and other 
magnetic confinement devices. The reactor design should provide for rapid 
replacement of damaged walls and electrodes. The desirable first wall neutron power 
flux would probably be about 2-5 MW m -2 . Higher values could cause heat removal 
and wall lifetime problems, and lower values would raise the cost of electricity. The 
anodes should be partially shielded from radiation, and the electrodes and walls 
must be actively cooled. Low-activation materials, such as purified vanadium- 
chromium alloys and SiC, should be used if feasible. 

One possible materials combination for an MEPC reactor is vanadium alloy 
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structure with a circulating Flibe (LiF-BeF 2 ) molten salt blanket and coolant. Use of 
the molten salt would avoid the reliability problems of high-pressure coolants and 
the fire hazard of lithium. A small amount of lithium could be used to cool first-wall 
hot spots, if needed. The irradiated vanadium alloy could be cooled for tens of 
years, then remelted to remove the volatile long-lived isotopes (mostly 42 Ar) and 
recycled [195]. With careful choice of materials, an MEPC reactor would not require 
active safety systems to prevent prompt off-site fatalities in the event of a serious 
accident. 


6. Conclusions and recommendations 

(i) Purely electrostatic plasma confinement with spatial and temporal focusing 
may be wroth re-examination (section 1.1.4). 

(ii) Due to magnetic reflection and electrostatic acceleration effects, the central 
plasma density in an MEPC device may be an order of magnitude higher than the 
electron density in the ring cusp region. Cold electrons from ionization of neutral 
gas tend to screen out the anode potential, but they are probably removed from ring 
cusps by diocotron oscillations (sections 2.5, 2.6). 

(iii) Although anomalous transport was observed in the ATOLL toroidal 
quadrupole device, the plasma loss rate in the Jupiter-2M linear ring cusp 
experiment was only about twice the classical rate, improving at higher density. The 
reason for the good stability of the Jupiter-2M plasma is not yet understood, but 
might be related to its method of electron injection, its better symmetry, or its large 
volume of uniform field-free plasma (section 3.7.3). 

(iv) Based on present experimental results, a linear set of ring cusps appears to 
be the most favourable configuration for an MEPC reactor. 

(v) Several decades of small, cheap, pulsed experiments have demonstrated the 
basic principles of MEPC, but these experiments have been marred by weak or 
pulsed magnetic fields, low voltages, weak plasma heating, incomplete diagnostics 
and inadequate computer modelling (section 3.7.4.). Further progress depends upon 
increasing R, B and <f> A in a steady-state device, so that n, T and r E can approach 
the thermonuclear regime. A new experiment is needed to demonstrate the 
attainment of interesting parameters (n > 3 X 10 1S> m -3 , T > 5 keV) in a steady-state 
device with intermediate size (R = 2 m, L ~ 10 m) and reasonable applied fields 
(B«5T, <£ A ~100kV). 

(vi) For an MEPC reactor with near-classical transport, values of Q = 10 are 
predicted at voltages 4> A ~ 400 kV (section 4.1). Q. could be enhanced by partial 
confinement of alpha particles. 

(vii) The Polywell tm and modified Penning trap concepts are just beginning to be 
studied. They need more research to determine the attainable degree of focusing 
and to explore related concepts that may lead to attractive fusion reactors (sections 
2 . 10 , 211 ). 

(viii) Radiation shields and voltage dividers should be used to prevent high- 
voltage breakdown. Voltage holding at 400 kV has been demonstrated in a magnetic 
mirror geometry by the PSP-2 rotating plasma experiment (section 5.3). 

(ix) MEPC devices may find applications as (section 4.2): electrical power plants; 
heavy ion sources; neutron sources for materials testing, medical therapy and 
isotope production. 
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(x) The main issues of concern for MEPC reactor development are: electron 
transport, plasma purity and electrode alignment and voltage holding. Future 
research should concentrate on resolving these issues. 
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